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EFFECT OF VELOCITY PROFILE DISTORTION IN CIRCULAR TRANSVERSE- 
FIELD ELECTROMAGNETIC FLOWMETERS 
by Norman C. Wenger 
Lewis Research Center 

SUMMARY 

The theory of circular transverse -field electromagnetic flowmeters is extended to 
include the effect of distortion from axial symmetry in the fluid velocity profile . A gen- 
eral formulation is developed for computing the flowmeter sensitivity for any given vel- 
ocity profile. It is shown how to best select the flowmeter pipe wall thickness and elec- 
trical conductivity to minimize error in the flow measurement due to uncertainties in the 
velocity profile. The ’’fully developed” velocity profile and corresponding flowmeter 
sensitivity are then computed for a Hartmann number range from 1 to 1000. Numerical 
results are presented in tabular form for a range of flowmeter parameters that is nor- 
mally found in liquid metal applications. These results enable the experimenter to cor- 
rect flow measurements for the effect of velocity profile distortion produced by the flow- 
meter itself. 


INTRODUCTION 

Electromagnetic flowmeters have been used for many years to measure the flow rates 
of a great variety of electrically conductive fluids including blood. The advent of nuclear 
reactors and associated energy conversion devices that use liquid metal coolants has 
caused renewed interest in electromagnetic flowmeters. 

Electromagnetic flowmeters use Faraday’s law as the basic principle of operation. 
The motion of a fluid throt^h a stationary magnetic field induces a potential gradient in 
the fluid. The value of this potential gradient, which is indicative of the flow rate, can be 
determined by measviring the potential difference between electrodes which are in elec- 
trical contact with the fluid. If the fluid has a high electrical conductivity, the potential 
gradient will result in a large induced cimrent in the fluid. The induced current produces 
a magnetic field which can also be measured and used to indicate the flow rate. Flow- 



meters which use this latter principle are generally referred to as induced-field flow- 
meters . 

Althot^h there are many t 3 ^es of electromagnetic flowmeters in use, the most com- 
mon tj^e for liquid metal applications and the only type that will be considered in this re- 
port is the transverse -field flowmeter shown in figime 1. The flow is contained in a cir- 
cular pipe which is normally made from nonmagnetic materials such as 316 stainless 
steel. A permanent magnet or direct current electromagnet is used to produce the uni- 


/-Known magnetic fiux density 



Figure 1. - Transverse-field electromagnetic flowmeter. 

form transverse static magnetic field. Electrodes for detecting the induced potential are 
located on a pipe diameter which is perpendicular to both the pipe axis and magnetic field. 
If the pipe wall is conductive, which is the most common situation found in practice, the 
electrodes are fastened to the outer wall of the pipe as shown. If the pipe wall is noncon- 
ductive, however, the electrodes must penetrate the wall to be in electrical contact with 
the fluid. 

There are several reasons for the popularity of this type of flowmeter. In addition to 
its basic simplicity, it does not require that the pipe be cut open during installation if the 
pipe wall is conductive since nothing is inserted directly into the flow stream. This is 
particularly advantageous in applications involving corrosive liquid metals where the 
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opening of a line to insert a flowmeter is both a possible source of contamination to the 
liquid metal and a hazard to the personnel involved. 

The main problem with this and all other types of electromagnetic flowmeters is 
determining the relation between the measured potential and the average fluid velocity 
or flow rate. Early studies of this problem by Williams (ref. 1), Thiirlemann (ref. 2), 
and Kolin (ref. 3) considered the case where the pipe wall was nonconductive and where 
the fluid velocity profile was symmetric about the axis of the pipe. It was found that the 
output potential is simply equal to the product of the pipe inner diameter, the applied 
transverse magnetic flux density, and the average fluid velocity. This resialt is valid for 
any velocity profile as long as the profile has axial symmetry. 

Several years later, Elrod and Fouse (ref. 4) considered the case where the pipe wall 
is conductive. They foimd that the output potential for this case is also proportional to the 
average fluid velocity as long as the velocity profile has axial symmetry. For the same 
flow rate, the outpvit potential is reduced from that foimd for a nonconducting pipe by a 
factor which involves the fluid and wall conductivity and the wall thickness. 

In many cases found in practice, the velocity profile of the fluid is not axially sym- 
metric. This asymmetry can be due to two effects. If^the flowmeter were located imme- 
diately downstream from an elbow or obstruction in the flow line, the velocity profile can 
be greatly distorted from axial symmetry. The obvious Solution of movii^ the flowmeter 
further downstream may not be possible in many cases. The second cause of velocity 
profile distortion from axial symmetry is the flowmeter itsej.f . If the fluid is highly con- 
ductive, as in the case of liquid metals, the large induced cvii'rent interacts with the ap- 
plied magnetic field and produces a force on the fluid. This force, which is nonsymmet- 
rical, is normally quite large and severely distorts the velocity profile from axial sym- 
metry. This distortion can be reduced by decreasing the intensity of the applied magnetic 
field. However, the output potential, which is approximately proportional to the applied 
magnetic field intensity, will normally vanish into the ambient electrical noise before the 
distortion becomes negligible. 

Considerable effort has been expended to determine how the velocity profile affects 
the output potential and, hence, the measurement of the average fluid velocity. A general 
formulation has been developed by Shercliff (refs . 5 and 6) showing the relation between 
the output potential and the velocity profile. This formulation is valid for any velocity 
profile, but it is limited to the case where the pipe wall is nonconductive. In the first part 
of this report, this formulation is extended to the case where the pipe wall is conductive. 
This formulation allows the experimenter who has some knowledge of the general shape of 
the velocity profile to compute a more accurate relation between the output potential and 
the average fluid velocity. It also shows how to best select the pipe wall thickness and 
conductivity to minimize the effect of uncertainties in the velocity profile. This formula- 
tion is most applicable when the fluid is a poor conductor so that the velocity profile in 
the flowmeter depends on conditions upstream from the flowmeter and may be predictable. 
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In the second part of this report, the effect of self-induced velocity profile distortion 
by the flowmeter itself is considered in detail. The exact solution for the velocity profile 
has been obtained by Ufland (ref. 7), Uhlenbusch and Fischer (ref. 8), and Gold (ref. 9) 
for the case where the pipe wall is nonconductive. The solutions are in the form of an in- 
finite series involving modified Bessel fimctions of the first kind. The significant figure 
requirement for summii^ these series exceeds that of modern computers, however, im- 
less the Hartmann nmnber M, which is the ratio of the electromagnetic forces to the vis- 
cous forces in the flxiid, is less than about 20. The range of Hartmann numbers encoim- 
tered in liquid metal applications is typically 50 to 500. 

An asymptotic solution for the velocity profile at large Hartmann numbers has been 
obtained by Singh and Noriboli (ref. 10) for the nonconductive pipe wall case. Although 
their results seem reasonable, their solution cannot be formally justified. The Bessel 
fxmctions in the series solution were replaced by an asymptotic form valid only when the 
argument, which is large but finite, is much larger than the order of the Bessel functions. 
The series was summed, however, from order zero to infinity (ref, 9). 

Some work has been done on the conductive pipe wall case when the wall has finite 
thickness, Shercliff (ref, 11) obtained a solution for the velocity profile and induced out- 
put potential usii^ an iterative technique. These solutions are in the form of power series 

2 

in the Hartmann number up to the term. M . The series give useful numerical values for 
the velocity profile and output potential only for M less than 2. 

Considerable work has been done on the conducting pipe wall case usit^ the "thin 
wall" approximation. In this approximation, the pipe wall is assumed to be sufficiently 
thin so that the effect of the wall can be incorporated into a boimdary condition at the pipe 
inner radius. Using this approximation, Diara, Tajima, and Matsushima (ref. 12) ob- 
tained a solution in infinite series form which also suffered from the same numerical dif- 
ficulties as did previous nonconductive wall solutions for Hartmann numbers greater than 
about 20. Also using the thin wall approximation, Shercliff (refs. 11 and 13), Chang and 
Lvmdgren (ref. 14), and Gold (ref. 15) obtained approximate asymptotic solutions at large 
M for the velocity profile and average fluid velocity. These solutions do not accoimt fully 
for the effect of the bovmdary layer at the pipe wall. Moreover, only Shercliff (refs. 11 
and 13) computed the output potential. The output potential solution is approximate and 
holds only over a limited range of Hartmann number and pipe wall conductivity. 

This report presents a solution for the output potential as a function of the average 
velocity for the fully developed velocity profile. The pipe wall is assumed to be conduc- 
tive and of arbitrary thickness . The solution is obtained using a variational principle for 
magnetohydrodynamic channel flow (refs, 16 and 17) and the Ritz approximation technique. 

The problem is also solved exactly using a procedure similar to that used in refer- 
ences 7 to 10 and 12. The exact solution, which is in infinite series form, also suffers 
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from numerical difficulties at large M. It is useful, however, at low M for testing the 
variational solution which is approximate. 

Sufficient numerical results from the variational solution are presented so that the 
experimenter, by interpolatir^ between the computed points, can obtain accvirate values 
for correcting meas\irements for the effect of velocity profile distortion due to the flow- 
meter itself. The parameter values covered are as follows: Hartmann number - 1, 2, 5, 
10, 20, 50, 100, 200, 500, and 1000; ratio of pipe wall to fluid conductivity - 0, .1, .5, 
1. , 1.5, and 2.; and ratio of outer -to -inner pipe diameter - 1., 1.05, 1.15, 1.25, and 
1.35. 

THE MODEL - BASIC EQUATIONS AND BOUNDARY CONDITIONS 
Preliminary Considerations 

The motion of a conducting fluid through an electromagnetic flowmeter is a complex 
process. The velocity profile of the fluid at the entrance to the flowmeter is determined 
by conditions upstream from the flowmeter. As the fluid approaches the magnetic field 
region, electric currents are induced in the fluid which react with the applied magnetic 
field and produce a force on the fluid. If the fluid is a poor conductor or the magnetic 
field is weak (i.e. , low Hartman number), this force is small and the fluid velocity pro- 
file remains essentially imchanged when passing through the flowmeter. However, if the 
Hartmann number is large, the velocity profile changes significantly with distance and 
time as the fluid travels through the magnetic field region. If the applied transverse 
magnetic field is uniform along the pipe axis and of sufficient extent, the velocity profile 
will eventually reach an equilibrium profile or "fully developed" profile. From this 
point on the profile does not change with distance or time until the fluid nears the down- 
stream end of the magnetic field region. After the fluid exits the flowmeter, the velocity 
profile eventually relaxes to axial symmetry barring any obstruction or elbow in the line . 

The ideal point to connect the electrodes for measuring the induced potential is after 
the velocity profile is fully developed. This choice has several advantages. If the elec- 
trodes were located at a point further upstream, the analysis becomes very difficult since 
the velocity profile at that point is dependent on upstream conditions which may be un- 
known. Furthermore, in this entrance region the induced potential varies with distance 
along the pipe. Consequently, the output potential depends not only on the velocity pro- 
file in the plane where the electrodes are attached but also on the velocity profile at 
points upstream and downstream as well. Any analysis of this region must, of necessity, 
be three-dimensional. In the fully developed region, however, there are no variations 
with respect to distance along the pipe so that a two-dimensional analysis is adequate. 
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Another factor that must be considered is whether the flow is turbulent or nonturbu- 
lent. Many, if not most, flows in liquid metal practice are turbulent. It has been foimd 
experimentally, however, that the imposition of a magnetic field transverse to the flow 
will suppress turbulence and raise the critical Reynolds number from about 2000 to a sub- 
stantially higher value. If the magnetic field strength is sufficiently high, the flow in the 
fully developed region can be made nonturbulent. This provides a great simplification in 
the analysis since the time varying terms in the fluid dynamic equations need not be con- 
sidered. 

Thus, to apply the resvilts of this analysis it is necessary to know the required dis- 
tance or entry length to achieve a fully developed profile for proper selection of the mag- 
net dimensions and electrode placement. And it is necessary to determine the critical 
Reynolds number to ensure nonturbulent flow. 

The required entry lei^h has been considered both theoretically and experimentally 
by Shercliff (refs. 6, 13, and 18) for the nonconductive pipe wall case. By usii^ a 
Rayleigh type of approximation, Shercliff fovind that for nonturbvilent flow the entry length 
L is of the order of 


L « Re -2- (1) 

M 

where Re is the Reynolds number, a the pipe inner radius, and M the Hartmann num- 

A 2 

ber which is assumed to be large. Using Re = 10 and M = 10 , which are typical val- 
ues found in practice, results in an entry length of 100 pipe radii. Experimental values 
for the entry length (ref. 13) are generally about one-half those given by equation (1). For 
the case of a conductive pipe wall the computed entry length is reduced by the factor 
l/yM(R - 1) (ref. 6) giving 


L « Re 2 (2) 

7M^(R - 1) 

where y is the ratio of the wall to fluid conductivity and R the ratio of the pipe outer - 
to-inner radii. Equation (2) is based on the assumption that 7M(R - 1) » 1. For the 
same values of Re and M used previously, the predicted entry length is 5 pipe radii for 
R = 1.2 and y = 1. An experimental investigation of the entry length has not been made 
for the conducting wall case nor has the problem of a turbulent entry flow been considered 
for either case. These shoiild be worthwhile areas to explore in future investigations. 

The required magnetic field strength necessary to suppress turbulence has received 
considerable attention, however. The experimental work of Shercliff (ref. 6) and 
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Loeffler, Maciulaitis, and Hoff (ref. 19) has indicated that the critical Reynolds number 
Recrit ^ circular pipe with a nonconductive wall is given by 

Hecrit = 250 M (3) 

for large M. No resrdts are available on the critical Reynolds number for the conductii^ 
pipe wall case. However, it is expected that the increased electric current in the fluid 
which results when the pipe wall is conductive will give greater stability to the flow and 
increase the critical Reynolds number from that given by equation (3) . Experimental ver- 
ification of this effect is also a worthwhile area for future investigations . 

Thus, to use the results based on the velocity profile in the fully developed region, it 
is necessary that the magnetic field be at least a few entry lengths in extent along the pipe 
axis and that the electrodes be attached at least an entry lei^th from the upstream end of 
the flowmeter. Furthermore, the magnetic field strei^th should be sufficiently large so 
that the Reynolds number is less than the critical Reynolds number to ensure nontiirbulent 
flow. 


The Model 

The model to be analyzed is shown in figure 2. The pipe has inner radius a and 
outer radius b and its axis coincides with the z -coordinate axis. The flow is in the 
positive z -direction. The magnetic permeability of the pipe wall is assumed to equal that 
of free-space. The applied magnetic field is uniform (independent of x, y, and z) and 
is alined parallel to the x-axis. If the pipe wall is conductive, the induced output poten- 
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tial is measured between points A and D. K the pipe wall is nonconductive, the poten- 
tial is meastmed between points B and C. In all cases it is assumed that the electrodes 
draw no current from the flowmeter and that the region beyond the pipe outer radius is an 
insulator . 


Basic Equations and Boundary Conditions 

The basic equations to be used are the standard magnetohydrodynamic equations for 
steady state, fully developed, incompressible, nonturbvilent flow. These equations, which 
consist of Maxwell’s equations, the continuity equation, the momentum transport equation, 
and the generalized Ohm's law, are given by 


V X E = 0 

(4a) 

V X B = PqJ 

(4b) 

V • B = 0 

(4c) 

V • E = ^ 

(4d) 



V . V = 0 

(4e) 

V)V = -Vp + J X B + 77V^V + PgE 

(4f) 

J = aj(E + V X B) + pgV 

(4g) 


where E, B, J, p^, €^, and are the electric field intensity, magnetic flux density, 
electric current density, electric charge density, electric permittivity of free -space, and 
magnetic permeability of free-space, respectively; and V, p^^, 77 , and p are the 
fluid velocity, density, viscosity, electrical conductivity, and pressiire, respectively. 

(All symbols are also defined in appendix A.) It has been assumed that the permittivity 
and permeability of the fluid are the same as that of free-space and that the fluid proper- 
ties n ,71, and (Jf are constant. 

It can be shown by using t 5 q»ical values for liquid metal properties that the terms in- 
volving the electric charge density p^ in equations (4f) and (4g) are negligible compared 
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with other terms in these equations (ref. 14). The term involving in equation (4d) 
must be retained, however. 

Some general properties of the solution of the basic eq\iations can be found without 
actually solving the equations . Since the pipe cross section and applied magnetic field do 
not vary along the z-axis, all quantities in the basic equations are independent of z ex- 
cept for the pressure. The pressure consists of the sum of two terms, one linear in z 
and the second involving an imknown function of x and y (ref. 20). Hence, the axial com- 
ponent of the pressure gradient is a constant. It can also be shown that the flmd velocity 

has only a z-component and the current density J has only x- and y-components. 

z 

Furthermore, the total magnetic field has only two components - the applied field in 

the X -direction and an induced field in the z -direction. 

Using these general properties of the solution, the basic equations can be reduced to 

a more manageable set. One approach is to combine all the basic eqioations yielding one 

fourth-order equation in V_. A second, and more preferable, approach which results in 

z 

two coupled second-order elliptic equations is to eliminate all dependent variables other 

than V and either the induced potential or the induced magnetic field, 
z 

First, consider retaining and the induced potential. The induced potential U 
can be defined as 


E = -VU (5a) 

Thus, equation (4a) is satisfied for any U. Next, taking the divergence of equation (4b) 
reveals that 

V . J = 0 (5b) 

Using this result along with equation (4g) with the p^V term omitted as discussed pre- 
viously then gives 


V^U = V • (V X B) (5c) 

The right side of equation (5c) can be simplified by remembering that V has only a 
z-component and that B has only x- and z-components. Thus, 


9^U , 9^U-p 


dx‘ 


9y 


9V. 


9y 


(5d) 


is the first of the coupled second-order equations involving only V^, and U. The second 
equation can be obtained by combining equations (4f), (4g), and (5a) to eliminate E and J 
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and then taking the z -component of the remaining vector equation givii^ 




ax‘ 


^^^z _ 1 3p °"f^o ^ 




3y" 


n 3 z V dy 


(5e) 


An equivalent set of coupled second-order equations involvii^ only V_ and B can- 

z z 

be obtained in an analogous manner. Takii^ the curl of equations (4b) and (4g) and com- 
bining gives 


V X V X B = 


fioCTf'ly X E + V 


X (V X B)] 


(5f) 


Expanding both sides of equation (5f) using standard vector identities and using equa- 
tions (4a) , (4c) , and (4e) and known properties of the solution to simplify the result give 


0^B 


ax" 


a^B 


ay 




(5g) 


The second equation can be obtained by combining equations (4b) and (4f) to eliminate J 
and then taking the z-component of the resviltii^ vector equation givii^ 


aV, 


ax" 


aV 


aB^ 
o z 


z _ ^ ^ 

gy2 77 az t 7 Pq ax 


(5h) 


The equations just derived apply, of course, only in the fluid. Equations which de- 
scribe the electromagnetic field in the pipe wall region can be obtained from the previous 
sets by setting the fluid velocity equal to zero. Thus, for the first set which involves only 

V and U the corresponding equation for the wall region is 
z 


a^u ^ a^u _ Q 

ax2 0y2 


(5i) 


Likewise, for the second set which involves only V„ and B_ the corresponding equation 
' z z 
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for the wall region is 


dx^ dy^ 


= 0 


(5j) 


In addition to these basic equations, appropriate boundary conditions must be speci- 
fied to determine the solution uniquely. These conditions are 


= 0 


r=a 


(6a) 


U = U (6b) 

r=a” r=a"^ 



a 

r=a w 


VB • £ 
z r 


r=a 


CTj VU 




lr=a 


lr=a 


= B 


lr=a" 

lr=a^ 

VU - a 1 

= 0 


r=b 
= 0 


r=b" 


(6c) 

(6d) 

(6e) 

(6f) 

(6g) 


where a^. is the vmit vector in the radial direction, is the electrical conductivity of 

the wall, and I _ refers to evaluating the quantity preceding it at a radial point just 
lr=a" 

slightly less than r = a. 

The boundary conditions require the following; 

(1) The fluid velocity must be zero along the pipe inner radius (eq. (6a)) . 

(2) The 6 or circumferential component of the electric field must be continuous 
across the fluid - pipe wall interface at r = a (eqs. (6b) and (6c)). 

(3) The radial component of the electric current must be continuous across the fluid - 
pipe wall interface at r = a (eqs . (6d) and (6e)) . 
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(4) The radial component of the electric current must vanish at the outer pipe radius 
r = b (eqs. (6f) and (6g)). 

Note that although there are only four stated boundary conditions they result in seven 
equations since equivalent forms involvir^ both U and are given. 

When solving the basic equations and associated boundary condition equations, it is 
convenient to reexpress the equations in the cylindrical coordinate system and to use di- 
mensionless coordinates and dimensionless variables . This can be achieved by intro- 
ducing the following quantities; 

Dimensionless potential; 


W = 


U 


B aV^ 
o o 


Dimensionless velocity; 


V = 



Hartmann number; 



Magnetic Reynolds number; 


R = u a^V 
m f^o f o 

Dimensionless axial pressure gradient; 

P = 

"" Wo 9z 


(7a) 


(7b) 


(7c) 


(7d) 


(7e) 
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Dimensionless induced magnetic field: 



Ratio of wall to fluid conductivity: 


r = 



Dimensionless radius: 


P 


r 

a 


Ratio of pipe outer-to-inner radii: 



a 


(7f) 


(7g) 


(7h) 


(7i) 


Characteristic velocity for the flow: 


(7j) 


Substituting equations (7a) to (7i) back into the basic equations and boundary conditions 
th.en yields the desired equations in dimensionless form. The equations are grouped so 
that the first set involves only the velocity and induced potential and the second set in- 
volves only the velocity and induced magnetic field: 


ap2 p 9p p2 g^2 9p p 00 


(8a) 


p3w^cos_e 

3p^ ^ de^ \ ^P P 8^ y 


p<l 


(8b) 


13 


(8c) 


Kp<R 


dp 


2 p 9p 2 gg2 


p=l‘ 


= 0 


W 


= w 


P=1‘ 


3=1 


+ 


(8d) 

(8e) 


9W 

dp 


= r 


p=l‘ 


9W 

dp 


p=l^ 


(8f) 


8W 

9p 


= 0 


3=R 


(8g) 


^ 1 3B ^ J_ 8^ e W . sin_e m 


8p 


2 P 9p n2 9p2 


m 


dp p dd 


(9a) 




8p2 P p^ de^ 


R 


m 


dp 


de, 


(9b) 


Kp<R 


dp 


2 P 3p p,2 g^2 


(9c) 


= 0 


p=l' 


(9d) 


r 


^1 

dp 


3=1 


9B | 

dp 


IP= 


(9e) 


B 


= B 


|p-l- 


P=1 


(9f) 
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Either the set of equations (8a) to (8g) or the set (9a) to (9g) completely specifies the 
problem. 


INDUCED POTENTIAL FOR AN ARBITRARY VELOCITY PROFILE 

In this section, a formulation is developed for computing the induced potential for a 
given but arbitrary velocity profile. The solution is obtained by using the dimensionless 
equations involving the velocity and induced potential (8a) to (8g). Equation (8b), which 
relates the fluid velocity to the forces on the fluid, need not be solved, however, since the 
velocity is assumed to be given. It is also assumed that the given velocity profile vanishes 
at the pipe wall so that equation (8d) is satisfied. 


Induced Potential 

The solution for the induced potential can be obtained in a straightforward manner by 
usii^ a Green's ftmction. The Green’s fimction G(p, ^ |pq, 0q) can be defined as the solu- 
tion of the following equations; 


a'^G ^ 1 dG ^ 1 
3p2 P ^P p2 

9 G _ _ 1 0 < p < R 

de^ Po 

(lOa) 


G _ = G ^ 

p=l p=l 

(10b) 


9G = 9G 

,+ 

p=l p=l 

(10c) 


^ =0 
9p 

p=R 

(lOd) 


where 6( ) is the Dirac delta function. 

By properly combining equations (8a) to (8g) with equations (10a) to (lOd) and using 



Green’s second identity it can easily be shown that the dimensionless potential W in the 
wall region is given by 




p,9) + 


cos 6 


o 3G 
90_ 




^(Po^VPo^^Po^^^o l^P^R ( 11 ) 

The potential in the fluid region 0 ^ p ^ 1 is also given by equation (11) with the 1/y 
factor omitted. 

The solution for the Green’s function is most easily obtained by using Fourier expan- 
sion techniques. Since equations (10a) to (lOd) reveal that G is an even function of 
9 - 9^, let 


OO 


G(p,e|po.®o) " Gn(P,Po)cos n{9 - 9^ 
n=0 


(12a) 


Substitutii^ equation (12a) into (10a) and usir^ the orthogonal properties of the functions 
cos n(0 - 9^) n = 0,1, 2, . . . over the interval of 9-9^ from 0 to 2ir show that 
Gn(p,Po) satisfies the equation 


where 


0% . 0G„ 

n ^ ^ n 


9p 


P 9p 



— - — 5(p - Pq) 

’^Po^n 


n 



if n = 0 
if n / O 


(12b) 


Since the solution of interest is the induced potential in the wall region 1 < p :< R 
due to the fluid motion within the pipe 0 < p ^ 1 it is only necessary to solve equation 
(12b) for 0 ^ p ^ 1 and l < p^ < R. The solution of equation (12b) subject to the bound- 
ary conditions (10b) to (lOd) is given by 


Go(p,Po) = 0 


(12c) 
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and 


Gn(p,Po> = 


rp 


n/ n 

• [Po 




n7r[(R^" + 1) + y(R^" - 1)] 


n^O 


0 ^ p 


1 SPo = R 


(12d) 


Thus, substituting equations (12c) and (12d) back into equation (12a) completes the solu- 
tion for the Green’s function. Combining the Green's function solution with equation (11) 
then yields a formula for computing the potential at any point in the pipe wall due to any 
fluid velocity profile. 

In practice, only the potential difference between points A and D is normally mea- 
sured (see fig. 2). This potential difference is given by 

Wad = W(R, jr/2) - W(R, -jr/2) (13a) 


Combining equations (11) and (13a) and using the Green’s function solution just obtained 
give 



4R" cos (n - l)e^ 
77 [(r 2“ + 1) + y(R^'^ - 1)] 


sin^V(p^,e^)p^dp^ 


d0. 


(13b) 


Flowmeter Sensitivity 

Rather than consider the absolute induced potential difference for a given flow condi- 
tion, it is more convenient to define a sensitivity for the flowmeter in terms of the poten- 
tial difference per imit of flow. It is also customary to include the applied magnetic field 
strength and the pipe inner diameter in the expression for the flowmeter sensitivity. 

Thus , the flowmeter sensitivity S can be defined as 

U 

S = — (14a) 

2a V„_ 


^If the pipe wall is an insulator, the same definition is used for S with and 

Wad ^®placed with and Wg^,, respectively. 
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where the mean or average flow velocity, is 


V. 


= — /* V (r, 0)r dr d0 
TTQ 2 Jq J 0 


zm 


Tra 


In terms of the dimensionless variables, S is given by 


S = 


w 


AD 


( 1 ) 


2V 


m 


(14b) 


(14c) 


where 


Vm=- /*y^’'v(p,0)pdpde (14d) 

“ 7T •'0 


Thus, the flowmeter sensitivity for any velocity profile is obtained by combing equa- 
tions (13b), (14c), and (I4d) giving 


S = 



(15a) 


where 


K(p, 0) 



n-1 

(R^*^ + 1) + y(R^" - 1) 


cos (n - 1)0 sin — 
2 


(15b) 


Equation (15a) is the desired result of this section. It shows that the flowmeter sensi- 
tivity is the ratio of the weighted average of the velocity over the pipe cross section to the 
average velocity with a weighting factor of K(p, 0) . 

Many of the early results on flowmeter sensitivity can be obtained quite simply from 
equations (15a) and (15b). First consider the case where the pipe wall is nonconductive. 
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This corresponds to setting y = = 0 and R = b/a = 1. The last step of setting R 

equal to 1 is a mathematically convenient way of moving the electrodes from points A 
and D to B and C without redoing the analysis (see fig. 2). It does not require the 
wall to be infinitely thin. If R = 1, y can actually assume any value since the term in- 
volving y in the expression for K(p,6) vanishes. Using these values in equation (15b) 
gives 

OO 

K(p,0) = ^ cos (n - 1)6 sin ^ (nonconductive pipe wall) (16) 

n=l 

This immediately shows that if the velocity profile has axial symmetry (i.e. , is indepen- 
dent of 6) then only the n = 1 term in K, which has unit magnitude, will contribute to 
the numerator in equation (15a). Thus, a flowmeter with a nonconductive pipe wall and an 
arbitrary but axially symmetric velocity profile has a sensitivity of 1 or 1 volt of signal, 
per meter per second of average flow, per weber per square meter of applied magnetic 
field intensity, per meter of pipe inner diameter. This result has been known for several 
decades . 

Next, consider the case where the velocity profile is still axially symmetric but 
where the pipe wall is conductive. Again, only the n = 1 term in the weighting factor 
K(p,6) contributes to the numerator in equation (15a). This time, however, the n = 1 
term has less than unit magnitude giving a reduced sensitivity of 

S= (17) 

(R^ + 1) + y(R^ - 1) 

This sensitivity is less than that for a nonconductive wall since some of the induced out- 
put potential is shunted by the conductive wall. This effect was first computed by Elrod 
and Fouse (ref. 4). 

The flowmeter sensitivity for a nonaxially symmetric profile cannot be computed 
until the profile is specified. However, much can be learned by examining the behavior 
of the weighting factor K(p, 6) over the pipe cross section. If K(p, 6) were a constant, 
the flowmeter would be ’’ideal; ” that is, its sensitivity would be completely independent 
of the velocity profile. 

Figures 3 and 4 show typical numerical values for the weightir^ factor K(p, 0) which 

O 

resvilt from summing the series in equation (15b). The values shown are normalized 

2 

The series expression for K(p,6) given by eq. (15b) can be expressed in closed 
form only for the case of a nonconductive pipe wall giving 

2 

K(p, 6) = - P (nonconductive pipe wall) 

1 + 2p^ cos 26 + p^ 

This result was also obtained by Shercliff (refs. 5 and 6). 


19 



Ratio of pipe wail to fluid conductivity, Ratio of pipe wall to fluid conductivity. 



Dimensionless radius, p 

(0 Cylindrical coordinate, 0 = 6rf>. Cyiindricai coordinate, 0 - 90<>. 

Figure 3. - Flowmeter weighting factor. Ratio of pipe outer-to-inner radii, R - 1.15. 
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Figure 4. - Flowmeter weighting factor. Ratio of pipe wall to fluid conductivity, y = 1- 00. 



with respect to the weighting factor at the center of the pipe to make the results more 
easily comparable. The value of the weighting factor at the center K(0,0) is 

2r/[(r2 + 1) + 7(r 2 - 1)] , the pipe shunting factor for axially symmetric flow. Because 

of symmetry considerations only one quadrant (0®:^ e 90^ need be considered. Fig- 
ures 3(a) to (d) show the normalized weightii^ factor as a fimction of radius at constant 6 
values of 0°, 30°, 60°, and 90°, respectively. The pipe outer -to -inner radii ratio R is 
1.15 and the wall to fluid conductivity ratio y is varied from 0 to 2. As shown, the 
normalized weighting factor is almost constant near the center of the pipe and drops off 
slowly near the wall for 6 less than about 60°. For 0 = 90°, however, the weighting 
factor increases quite rapidly with radius and, in fact, becomes infinite at the wall for 
y = 0 . The flowmeter output is , therefore , most sensitive to the flow at points near the 
wall at 0 « 90° and least sensitive to the flow near the wall at 0 « 0°. The large ratios 
of wall to fluid conductivity make the weighting factor closer to that of an ideal flow- 
meter which has a constant weighting factor. 

Figures 4(a) to (d) show additional numerical results for K(p, 0) only this time the 
ratio of wall to fluid conductivity y is held constant at 1 and the pipe outer -to-inner 
radii ratio R is varied. The results show that as the wall thickness increases the 
weighting factor becomes more nearly a constant. 

It can be concluded, therefore, that flowmeters with the highest ratio of pipe outer- 
to-inner radii and wall to fluid conductivity are the least sensitive to the shape of the 
fluid velocity profile . This insensitivity to the velocity profile is achieved only at the 
price of reducing the overall sensitivity since a thick, highly conductive pipe wall gives 
the largest shunting effect. Thus, the experimenter who can select the pipe wall thick- 
ness and conductivity must consider at least two effects in minimizing the flow measvire- 
ment error. The pipe wall must be svifficiently thick and conductive to minimize errors 
due to the uncertainty in the exact velocity profile, and at the same time the wall must 
permit a sufficiently large output potential to minimize the measurement error due to the 
ambient electrical noise . 

' As mentioned in the INTRODUCTION, the principal use for the preceding analysis is 
in the case of low Hartmann number where the velocity profile in the flowmeter is deter- 
mined by upstream conditions. In the following section, the case where the velocity pro- 
file is determined by the flowmeter itself is considered. 


INDUCED POTENTIAL FOR FULLY DEVELOPED VELOCITY PROFILE 

In this section, the entire set of basic equations will be solved simviltaneously to ob- 
tain the self-consistent solution for the fully developed velocity profile and the corre- 
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spending induced potential. From these results the flowmeter sensitivity can be easily 
computed. 

Two approaches will be followed. First, an exact series representation of the solu- 
tion will be obtained, and second, an approximate solution using a variational principle 
and the Ritz technique will be computed. The exact series solution can be summed 
numerically usii^ modern computers only for Hartmann numbers M less than 10 to 20, 
whereas the solution foimd using the variational principle can be evaluated at all Hartmann 
numbers. The exact series solution is of value, however, since it gives some insight on 
the accuracy of approximate solution in the lower Hartmann number range. In addition, 
it is much faster computationally to use the series solution where it applies. Finally, by 
using transformation techniques such as the Watson transform it may be possible to 
transform the exact series solution which converges slowly at high Hartmann nmnbers 
into one which converges rapidly. 


Exact Series Solution 

Sol ution of basic equations . - In obtaining the series representation of the solution it 
is convenient to work with equations (9a) to (9g) which involve only the velocity and induced 
magnetic field. Equations (9a) and (9b)' can be uncoupled by introducing the auxiliary fimc- 
tions f and g where 


V(p, e) = i 
2 


M 


M 


p cos 0 — p cos 6 

f(p, 0) + e ^ g(p, e) 


(18a) 


B(p, e) = 


R 

m 

2M 


M M ■ 

— p COS 6 — p COS 9 

2 f(p,9)-e^ g(p,9) 


P R 
o m 




p cos 6 (18b) 


Substituting equations (18a) and (18b) back into equations (9a) and (9b) reveals that f 
and g must satisfy the following equations: 


9p2 P p2 302 4 


(19a) 
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(19b) 

3p^ P ^p de^ ^ 

The general solution for f and g can be expressed as a Fourier series in 6. Sym- 
metry considerations show that both V and B, and, hence, f and g, must be even fimc - 
tions of 6 about 6 = 0. Thus, only cos n0 type terms need appear in the series. The 
coefficients in the Foimier series are, in general, linear combinations of the modified 
Bessel functions with argument Mp/2. Since modified Bessel fxmctions of the second 
kind diverge at zero argviment, they must be excluded. Thus, the solutions to equations 
(19a) and (19b) are given by 


OO 

«P.0)=^Vn(f 

n=0 


p) cos n6 


OO 

g(p,«) = y 

n=0 


M 


p) cos n0 


(20a) 


(20b) 


where and B^^ are constants to be determined. Combining the solutions for f and g 
with equations (18a) and (18b) give the general solution for V and B in the fluid region 

0 ^ pS 1: 



p cos 6 


'nl^nfyp) 


0 < p < 1 


(21a) 


r»-n 


B(p,0) 




cos ne 


m 2 


p cos 6 0 ^ p £ 1 


(21b) 
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In the pipe wall region 1 < p < R the induced magnetic field satisfies equation (9c) 
which is Laplace's equation. The solution of this equation which satisfies boundary con- 
dition (9g) can also be found using Fourier expansion techniques and is given by 


B(p,0) 


oo 

y C„(p" - r 2 V”) cos ne 
n=l 


1 < p < R 


(21c) 


where the are constants to be determined. 

The unknown constants A^, and Cj^ can be determined by requiring the general 
solutions given by equations (21a) , (21b) , and (21c) to satisfy the three remaining bo\md- 
ary conditions (9d), (9e), and (9f). Substitutii^ equations (21a) to (21c) into eqviations 
(9d) to (9f) and using the orthc^onal properties of the cos n0(n = 0, 1, 2, . . .) fxmctions 
over the interval 0 ^ 6 ^ 2ir and the integral definition of the modified Bessel fimctions 



give the following relations: 

= 0 n even (23a) 


n 


r 

n(R^" + 1) 


R P 
m o 






R P 


m o 


2M 



•? Z 

k— 1 ^3^5^ . • • 


(1 - R^*^)S 



n odd 


(23b) 


25 




26 



where 


Cl if i=j 

Lo if ij^i 


and the prime (') denotes differentiation with respect to the argument of the modified 
Bessel functions. 

These equations show that the (n = 1 , 3, 5, , . , which are the coefficients in the 
expansion for the induced magnetic field in the pipe wall, satisfy an infinite set of linear 
algebraic equations. The solution for the can be obtained by truncating the infinite 
set at a siafficiently large n and then solving the remainii^ finite set usii^ standard tech- 
niques. It should be noted that in this set of equations both the inhomogeneous terms, 
Tn(M/2), and the coefficients of the terms, S^(M/2), involve infinite series (see 
eqs. (23e) and (23f)). These infinite series, however, involve only the Hartmann number 
M and are independent of the conductivity and radius ratios y and R. Thus, the series 
need be summed only once for each M of interest. 

Once the have been determined, it is then quite straightforward to compute the 
using equation (23c) and the from eqtiation (23d). The velocity and induced mag- 
netic field can then be found from equations (21a) to (21c). Before presentir^ numerical 
results, expressions for the average velocity, induced potential, and flowmeter sensi- 
tivity will be obtained. In addition, the limiting case of low M will be considered. 

Average velocity and induced potential. - The flowmeter sensitivity, in terms of the 
dimensionless variables, is the ratio of the potential difference to twice the mean 

velocity (see eq. (14c)). Thus, to compute the sensitivity for the series solution 
just obtained it is necessary to relate the induced potential to the velocity and induced 
magnetic field. This is most easily done by first computing the electric field and then 
evaluating the line integral of the electric field from one electrode to the other to obtain 
the potential difference. Since the answer is independent of the path of integration, it is 
most convenient to use the circumferential path along the outer surface of the pipe wall. 

In the pipe wall region. Maxwell’s equations and Ohm’s law can be combined to give 

V X B = (24a) 

Since the path of integration corresponds to a constant radius, only 8U/30 need be eval- 
uated: 

= ^(V X B) • a. = — (24b) 

9 ^ Po^w 
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Reexpressing equation (24b) in terms of the dimensionless variables then yields 


^ = (24c) 

90 9p 

The potential difference can be computed by integrating equation (24c) with respect 

to 9 from -n/2 to tt/ 2 along the curve p = R: 


W 


AD 



de 


(24d) 


Combining this expression for with the general solution for B in the wall region 

(eq. (21c)) then gives 


W 


AD 


yR 


m 


oo 

z 


C„R" sin ^ 
2 


n=l,3, 5, . . 


(25) 


The average dimensionless velocity can be obtained by combinii^ the definition 

for the average (see eq. (14d)) with the general solution for V given by equation (21a) 
and performing the integration with respect to p and 9. This gives 


m 


oo 

E 

n=0 


(-1)"A 


n 


I (M ^ + 4n^) j2 IM\ 


n ' 





(26) 


where equation (23d) was used to eliminate the terms. The flowmeter sensitivity is 

then simply the ratio of (eq. (25)) to twice (eq. (26)). 

Asymptotic solution for small M . - Results for the limiting case of low Hartmann 
number can be obtained by expressing T^(M/2) and S^(M/2) in a power series in M by 
using the asymptotic form for the modified Bessel functions for small argument. Substi- 
tutii^ these expansions into the set of equations for the (eq. (23b)) then gives 
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(27g) 


Substituting these residts for the and into the general expressions for the in- 
duced potential difference W^j^, the average velocity V^, and the sensitivity S then 
gives 
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r 


w 




AD 


2 [(r 2 + 1) + y(R^ 



_ y(R^ - 1)M^ 

8 [(R^ + 1) + y(R^ - 1)] 



r o o 1 1 

R^M^ 

L(r2 + 1) + y(R^ - 1)J 

576 

T . . • ► 

[(R® + 1) + y(R® - 1)J ^ 


(28a) 


y(R^ - 1)M^ ^ 


m g 

48 

8[(r 2 + 1) +y(R^ - 1)] 




r 



2R 


^ _ mV 

(R^ + 1) + y(R^ - 1) 

+ > 

[(r 2 + 1) + y(R^ - 1)] 1 

1^ 576 

(R® + 1) + y(R® - 1)_ 

T . . . r 


(28b) 


(28c) 


This same result has also been obtained by Sher cliff (ref. 11) who used an iterative 
procedure. The first term in the expression for the sensitivity is recognized as the pipe 
wall shunting factor. The results show that the fully developed velocity profile at low 
M causes a decrease in the sensitivity. Equation (28c) gives useful numerical re- 
sults (within 1 percent) for M 4. Extending the results in the same manner including 
terms that vary as expands the useful range to only M :< 4. 5. In the next section, a 
numerical evaluation of the exact series solution is presented. 

Numerical results for exact s eries solution . - The infinite series for T^(M/2) and 
Sn|^(M/2) were summed usir^ an IBM 7094 computer and double precision arithmetic. 

The modified Bessel fxmctions were evaluated to at least 14 place accuracy. The prin- 
ciple difficulty in summir^ these series for Tj^(M/2) and Sjj^(M/2) is that the terms al- 
ternate in sign and the values of the sums are normally many orders of magnitude less 
than the magnitude of the largest term in the series. To obtain a minimum of 3 to 4 place 
accuracy for T^(M/2) and S^(M/2) the values for n, k, and M were restricted to 
where the values for T^(M/2) and S^(M/2) were at least 10“^ to 10"^® times the largest 
term in their respective series. This restricted n and k to be less than 19 for M of 
5 and 10, and less than 11 for M of 20. 

The limitation on the maximum values permitted for n and k for a given M imme- 
diately determines the maximum number at which the infinite set of equations for the 
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must be trimcated. A detailed study of these equations reveals that they are diagonally 
dominant, at least for M 20. The computed value of Cj, for example, was found to be 
almost independent of the number of equations used. Thus, the values obtained for the 

are accurate, but the number of that can be computed is limited. 

The next difficulty arose in evaluating the given by equation (23c) because of the 
infinite series. Often this series did not converge usii^ the finite number of avail- 
able. ""he nonconductive pipe wall case (R = 1) presented no trouble at this point since the 
series vanishes so that A^^ could be computed for any n. However, as the pipe wall 
thickness and conductivity and the Hartmann number are increased, the maximum number 
of that can be computed decreases since there are insufficient numbers of avail- 
able to sum the series. For example, usii^ the thickest wall (R = 1.35) and largest con- 
ductivity ratio of interest (y = 2.0) permits the A^ up to Aj^g to be calculated to 4 place 
accuracy at M = 10 but only up to Ag to the same accuracy at M = 20. 

Additional limitations occvirred in evaluating the velocity profile, the induced poten- 
tial difference W^jj, the average velocity V^, and the ratio W^/2V^ which is the 
flowmeter sensitivity. There were insxifficient numbers of A^ available to evaluate the 
series solution for the velocity profile (eq. (21a)) at M = 20 except for the case of a non- 
conductive pipe wall. As the pipe wall thickness and conductivity increased, the values 
for M had to be restricted in order to sum the series to 4 place acciaracy with the finite 
number of terms available. For thickest and most highly conductive pipe wall case con- 
sidered (R = 1. 35, y = 2.0), the velocity profile series covild only be summed for M ^ 12. 

The infinite series for the average velocity over the pipe cross section (eq. (26)) con- 
verged faster than the series for the velocity at any given point in the cross section. The 
series for could be evaluated to at least 5 place accuracy for M ^ 10 and to a max- 
imum error of a few parts in the third place at M = 20. The series for the induced poten- 
tial difference converged rapidly for M < 20 and presented no difficulty. 

The numerical results are presented in table I for M values of 1, 2, 5, 10, and 20. 
The values are accurate to the 3 places shown except at M = 20 where the third place 
may be in error. The values used for the ratio of the outer-to-inner pipe radii R are 
1.00, 1.05, 1.15, 1.25, and 1.35 and the values used for the wall to fluid conductivity 
ratio y are 0, 0.1, 0.5, 1.0, 1.5, and 2.0. The motivation for selectii^ these particu- 
lar values is presented in appendix B. 

Results are given for W^/P^, S(=W^/2Vjj^), and the ’’distortion factor” 

Fp . The distortion factor Fjj is related to the flowmeter sensitivity by 

S = 2R Fj 3(R, y , M) (29) 

(R^ + 1) + y(R^ - 1) 
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It is introduced to clearly separate the two causes for reduced sensitivity. The first is 
simply the electrical shunting effect of the induced output potential by the conductive pipe 
wall for the axially symmetric profile. This is the first term in equation (29), and it is 
independent of M. The factor accoimts for the change in sensitivity due to the dis- 
tortion of the velocity profile from axial symmetry by the applied static magnetic field. 

The results in table I show that the distortion factor which is 1 at M = 0, de- 
creases as M increases. This result is in agreement with the asymptotic solutions for 
small M given in the previous section. What is surprising, however, is that Fjj de- 
creases more rapidly with increasing M if R and y are small. This is not expected 
since small values for R and y produce the least distortion in the velocity profile at 
low M and, hence, are expected to yield the smallest change in F^. However, the pre- 
vious section, which considered the effect of an arbitrary velocity profile, showed that 
flowmeters with small R and y are the most sensitive to distortions in the velocity pro- 
file from axial symmetry. Thus, even though flowmeters with small R and y produce 
the least distortion for M < 20, they produce the largest change in the distortion factor 

^D* 

The remaining entries in table I show that both and decrease with 

increasing M and that the decrease is more rapid for large R and y. This effect is 
due to the increase in P^ which is proportional to the pressure drop per unit length 
along the pipe. As M increases, the velocity profile distorts from its Poiseuille profile 
resulting in a larger viscous drag force. In addition, if the pipe wall is conductive, some 
of the induced cxirrent, which must form a closed path, can return through the pipe wall 
rather than only through the fluid, thus giving a net electromagnetic drag force on the 
fluid. The increase in these forces must be balanced by a corresponding increase in P^ 
to maintain the same flow. The values for V^^/P^ can be used to predict the pressure 
drop across a flowmeter. However, it must be remembered that the tabulated values are 
only for the fully developed region. Appropriate amounts must be included for the entry 
and exit regions to obtain the total pressure drop. 

In the next section, an approximate solution is obtained using a variational principle 
and the Ritz technique . 


Approximate Solution Using Variational Principle and Ritz Technique 

Variational principl e. - An alternative approach for obtaining a solution is to convert 
the problem of solving the governing differential equations and associated boundary condi- 
tions into an equivalent variational problem. It can be shown (refs. 16 and 17) that the 
Euler -Lagrange equations and natural boiandary conditions of the functional F are the 
basic equations and boundary conditions (8a) to (8g) where F is given by 
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ew 

00 


>p dp d0 - yM 



p dp d0 



The functional F has the property that its first variation is zero. Thus, F is sta- 
tionary; that is, first-order changes in V and W about their true values produce only 
second -order changes in F. It can also be shown that of all functions that are continuous 
with piecewise continuous first derivatives, the particular functions for V and W 
that make F stationary satisfy equations (8a) to (8g) and, hence, are the desired solu- 
tions . 

A very important property of the functional F (refs. 16 and 17) is that its stationary 
value Fg|. is given by 



Since the dimensionless pressure gradient is a constant, the stationary value for F 
is proportional to the average fluid velocity in the pipe. This is a very useful result since 
the average velocity is one of the quantities needed to obtain the flowmeter sensitivity. 

The stationary value of F can normally be computed to good accuracy by usii^ approxi- 
mate fimctions for V and W. 

A careful study of equation (30) reveals that the stationary value of F is not an ex- 
tremum. However, if the class of admissible fimctions for V is restricted so that either 
V or dV /dp is specified at p = 1”, then it can be shown that the stationary value for F 
is a maximmn (refs. 16 and 17). Thus, if F is maximized over a subset of the entire 
class of admissible functions for V and W subject to this restriction, the resulting 
stationary value when used in conjunction with equation (31) will give a computed average 
velocity which is less than or equal to the actual average velocity. The proper restric- 
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tion to impose is to reqvdre all admissible fimctions for V to be zero at p = 1” which 
is boimdary condition (8d) . 

Ritz technique . - Approximate solutions for the velocity and potential will be ob- 
tained using the Ritz technique. In this technique, the velocity and potential are ex- 
pressed in terms of known functions of p and B that approximate the true solution but 
contain adjustable parameters ^ , . . . , Substituting this approximate solution into 

equation (30) and performing the integrations over p and 9 then leave F as a function 
of the adjustable parameters ^ , . . . , and the characteristic parameters of the 
model P^, M, y, and R. Assuming that the approximate solution for V vanishes at 
p = 1” for all values of . . . , the stationary value of F can be found by maxi- 
mizing F with respect to The corresponding values for . . . , at 

the stationary point, when substituted into the approximate fvmctions for V and W, will 
yield the closest apprcecimations to the velocity and potential that are possible for the 
class of functions used. 

It is convenient to express the general solutions for V and W by the followii^ 
Fourier expansions: 


N 

V(p, ^) = Yj 

n=0 

N 

W(p, ^) = 2 
n=0 


where the radial fimctions f^^ and g^^ are to be determined. Symmetry considerations 
show that the cos 2n0 and sin(2n + l)e functions with n = 0, .1, 2, ... °° are com- 
plete for V and W, respectively. Substituting expansions (32a) and (32b) into equa- 
tion (30) and performii^ the integration over B then give 


N 



34 


I 



Each of the radial functions f^^ and gj^ can be expressed in series form as 


^ ^niP*""' " = , N; 0 s p ^ 1 


(34a) 


i=l 


n 

'n<P> ' 


=n,p'"* 


n = 0,1,2, . . . , N; 0 ^ p s R 


(34b) 


3=1 


Since the functions p”(n = 0,1,2, . . . , «>) form a complete set with respect to functions 
that are continuous with piecewise continuous derivatives, the exact solution can be ob- 
tained by setting = i - 1, = j - 1, N = «>, = «>, and = <», by substituting 

equations (34a) and (34b) into equation (33) , and by finding the stationary point of F by 
solving 


.0 

n = 0,1,2, . . 

., N; 

i — 1,2,3, . . 

(35a) 

=0 

n = 0,1, 2, . . 

., N; 

] = 1,2,3, . . ., 

(35b) 


for and . The solution of this doubly infinite set of equations yields the exact 
solution to the problem. In practice, these equations are solved by first trxmcating the 
sets at some arbitrary upper limit and solving the resulting equations. The upper limit 
is then increased and the process is repeated tintil the sequence of solutions converges. 
Convergence to the correct solution is assured since the ftmctions used in the expansions 
form complete sets. 

The main difficulty with this technique is that, for large M, the number of eqviations 
to be solved of the form (35a) and (35b) become prohibitively large before convergence is 
achieved. The reason for this difficulty at high Hartmann number is that the functions 
fj^, for example, are normally very small for p < 1 - (1/M) (except for fQ) and become 
large and oscillatory with p for 1 - (1/M) < p < 1. It takes a large nimiber of functions 
of the form p*^ (n = 0, 1, 2, . . .) to represent a function with this behavior. 

It was found, however, that if the exponents and /3j^. in the trial functions (34a) 
and (34b) are also considered as adjustable parameters rather than fixed integers then 
convergence is achieved with a manageable number of equations. Thus, in addition to 
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solving equations (35a) and (35b), equations of the form 


3F 


3 O' 


= 0 


ni 


0F 


3^. 


= 0 


n] 


n = 0,1,2, 
n = 0,1,2, 


. , N; j = 1, 2, 3, . . • j Jjj 


(35c) 


(35d) 


must also be solved. Since the exponents are variables, completeness of the set of ex- 
pansion functions is no longer assured. However, this is not a serious drawback since 
much is known about the general form of the solution and various tests are available for 
checking the results. 

In order to further minimize the number of equations to be solved, several of the 
known properties of the solution were built into the trial functions. In the pipe wall re- 
gion, the radial functions g^^ which satisfy equations (8c), (8e), and (8g) are given by 


. , , p2n+l ^ ,j4n+2 p-(2n+l) ^ 

VP> = Zri; 2v° 


4„.2 ^ n] 

1 + j=l' 


n = 0,1,2, . . ., N; 1 < p < R 


(36a) 


Thus, the expression for gj^ given by equation (34b) is used only over the interval 
0 s p s 1 and equation (36a) is used for g^^ for 1 ;< p < R. In addition, the boundary 
condition which requires V to vanish at p = 1 " can be built into the trial functions quite 
simply by setting fj^(l) = 0 or 


=0 n = 0,1,2, ..., N (36b) 

i=l 

Equation (36b) ensures that the stationary point of F is a maximum. Thus, rather than 
solving equations (35a) to (35d) to obtain the stationary point, an alternative approach is 
to maximize F with respect to the and j3^. subject to the constraint 

given by equation (36b). 

Evaluation of solution . - In most applications of the Ritz technique, the adjustable 
parameters are restricted to the linear coefficients such as the A^^^^ and in equa- 
tions (34a) and (34b). The equations to be solved for the stationary point ((35a) and (35b)) 
are, therefore, linear in the adjustable parameters since f^^ and g^ appear in F only 
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to the second power. The solution is thus reduced to solving a system of sim\iltaneous 
linear algebraic equations. 

Allowing the exponents in equations (34a) and (34b) to vary still results in linear 
equations for the and Cj^. ((35a) and (35b)) but they are nonlinear in and 
In addition, the other equations that are introduced ((35c) and (35d)) are nonlinear in all 
variables . 

Rather than attempt to solve these nonlinear eqviations, it was found simpler to max- 
imize F directly. This was done by using the search technique developed by Powell 
(ref. 21) which finds the minimum of a function of several variables without calculating 
derivatives. (Maximizing F is the same as minimizing -F.) This technique was used 
primarily to determine how many adjustable parameters were needed for convergence as 
a function of M, R, and y. The problem was repeatedly solved for increasing values of 
N, I^, and where N is the upper limit in the summation over the angular index n 

(see eqs. (32), (33), and (34)) and 1^^ and are the upper limits in the summations 
over the radial indices (see eqs. (34a) and (34b)). The process was terminated when an 
increase in N, I^, and did not produce a corresponding increase in the maximum 
value for F. 

It was found that at a given Hartmann number, the largest number of parameters re- 
quired for convergence occurred when the pipe wall was nonconductive (i.e. , R = 1.0, 
y = 0) . The number of parameters was then fixed at the value required for convergence in 
the nonconductive pipe wall case when evaluatii^ results for other values of R and y at 
the same M. Table II shows the number of parameters used at each value of M. 

Although this technique for maximizing F worked well even if the initial guess was 
poor, it was too slow for computing all the results. Consequently, once the solution for 
the nonconductive pipe wall case at each M was determined, a modification of another 
technique developed by Powell (ref. 22) for solving systems of nonlinear algebraic equa- 
tions was used. This technique, which is a combination of Newton's method and the me- 
thod of steepest descents, was used to solve equations (35a) to (35d) as follows: 

(1) Using the best guess available for the exponents and equations (35a), 

(35b), and (36b) were solved exactly for A^^j^ and This could be done quite rapidly 

since these equations are linear algebraic equations in A^j and 

(2) The solution for the A^^ and with the corresponding values for and 

^nj used to evaluate 0F/9 q!^j^ and 9F/9j8jjj (left sides of eqs. (35c) and (35d)). 

(3) Powell's technique was then used to compute values for 01 ^^^ and that satis- 
fied equations (35c) and (35d). Each time new values of 9F/9(¥j^j and 9F/9j3^j were 
required steps (1) and (2) were repeated. 

(4) The iterative procedure was terminated when the sum of the squares of 9F/9(v . 

and less than a cutoff value which Avas typically to 10” 

(5) After completing steps (1) to (4) for given values of R and y, the process was 
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repeated for new values for R and y of interest at the same M. The new guess used 
in step (1) was the solution for the previous R and y. 

Once the stationary point for F and the values of the adjustable parameters at the 
stationary point were computed, the approximate solutions for V and W could be de- 
termined by combining equations (32a), (32b), (34a), (34b), and (36a), 

Numerical results . - Numerical results obtained using the variational principle and 
the Ritz technique are presented in table in. Values are given for the average velocity 
Vm/^o’ induced potential difference flowmeter sensitivity S, and distor- 

tion factor for M values (1, 2, 5, 10, 20, 50, 100, 200, 500, and 1000), R values 
(1.00, 1.05, 1,15, 1.25, and 1.35), and y values (0, 0.1, 0,5, 1.0, 1.5, and 2.0). The 
results for the distortion factor are also presented in figures 5(a) to (d) . 

At low values of M, the distortion factor decreases with increasing M, The rate of 
decrease is greatest for pipes with a thin, poorly conducting wall and least for pipes with 
a thick, highly conductive wall. These results are in exact agreement with those fovind 
previously from the exact series solution. As M continues to increase, however, the 
distortion factor reaches a minimum if the pipe wall is conductive, and then it begins to 
increase and to asymptotically approach 1 as M becomes infinite. If the pipe wall is 
nonconductive, the distortion factor continues to decrease with increasing M and 
approaches approximately 0. 924 as M becomes infinite. 

The results in figure 5 show that it is important to include the effects of velocity pro- 
file distortion by the flowmeter itself when makii^ precision flow measurements. If these 
effects are not considered, the indicated flow measurement will be less than the actual 
flow. If the pipe wall is nonconductive, this error increases with increasing M and ap- 
proaches (1 - 0.924/0. 9 24)xi 00 or 8.2 percent as M becomes infinite. If the pipe wall 
is conductive, the error reaches a maximum, typically 2 to 5 percent, in the Hartmann 
number range 10 to 50 and then decreases to zero as M becomes infinite. 

The reason for the behavior of Fjj with increasing M can- be seen by examining the 
velocity profiles. Figures 6 to 8 show the variation of V with p at fixed 6 values of 
0°, 30°, 60°, and 90°. The velocity is normalized with respect to V^, the dimension- 
less velocity at the center of the pipe. Three cases are shown - the nonconducting pipe 
wall (R = 1.00, y = 0), a thin, poorly conducting pipe wall (R = 1,05, y = 0,1), and a 
thick, highly conducting pipe wall (R = 1.35, y = 2.0). 

At M = 0, the velocity profile is axially symmetric and assumes the normal 
Poiseuille shape for all cases. As M increases, the velocity profile becomes flattened 
particularly along 6 = 0°. As M continues to increase, the profile also flattens along 
6 = 90° if the pipe wall is conductive, and the profile eventually approaches slug flow as M 
becomes infinite. If the pipe wall is nonconductive, however, the flattenii^ effect is not 
as pronounced along 6 = 90°, and the profile does not approach slv^ flow as M becomes 
infinite , 
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(c) Cylindrical coordinate, 0 = 60° (d) Cylindrical coordinate, 0 ^ 90°. 

Figure 6. - Velocity profile. Ratio of pipe outer-to-inner radii, R = 1.00; ratio of pipe wall to fluid conductivity, y = 0. 
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Dimensionless radius, p 

(c) Cylindrical coordinate. 0 = 60°. (d) Cylindrical coordinate, 9 = 90°. 

Figure 7. - Velocity profile. Ratio of pipe outer-to-inner radii, R = L 05j ratio of pipe wall to fluid conductivity, y = a 1. 




The reason for this flattening phenomenon is as follows . If the pipe wall is noncon- 
ductive, all the induced cvurrent must remain in the fluid since the current must form a 
closed path and there can be no current in the pipe wall. The net electromagnetic J x B 
force on the fluid is, thus, zero. Since J x B is not identically zero everywhere, the 
effect of this force is to distort the profile . The magnitude of this force is greatest near 
the center of the pipe where it opposes the flow and along the pipe wall near 6 = 0° and 
180° where it aids the flow. This force is minimal near the pipe wall at 6 = 90° and 270° 
since J and B are essentially parallel in this region. 

If the pipe wall is conductive, a fraction of the induced current in the fluid returns 
through the pipe wall. It was found that this fraction increases as M increases so that 
at high M virtually all the return current is in the pipe wall. In addition, it was found 
that the current density in the fluid becomes uniform over the cross section at high M 
giving a uniform electromagnetic force. The net effect is that the uniform electromag- 
netic and pressure forces balance each other in the core region resulting in a constant 
core velocity with negligible viscous forces. The viscous forces are only important in 
the thin bovmdary layer region where the velocity drops sharply from its core value to 
zero at the wall. 

The reason for the variation of the distortion factor with M is now clear. If M is 
small, the distortion factor decreases with increasing M since the velocity profile be- 
comes distorted because of the nonsymmetric electromagnetic forces. The distortion 
factor decreases rather than increases since the flow is beii^ shifted from a region where 
the sensitivity weighting factor is high (along 6 = 90°, 270°) to where it is low (along 
6 = 0°, 180°) (see figs. 3 and 4). As M continues to increase, the velocity profile goes 
to slug flow if the pipe wall is conductive, and thus the profile returns to axial symmetry. 
This causes the distortion factor to return to 1. If the pipe wall is nonconductive, the 
profile remains distorted from axial symmetry as M increases. The distortion factor 
for this case approaches 0.924 as M becomes infinite. 

Accuracy of results . - The last main consideration is to determine the accuracy of 
the numerical results just presented. In most applications of the Ritz technique, the num- 
ber of adjustable parameters is increased until the stationary value of the functional con- 
verges to a definite value. If the trial fimctions are expressed using a complete set of 
functions, convergence to the correct result is assvired barrii^ any error that may be in- 
troduced in the computation process. It is often difficult, however, to determine if con- 
vergence has actually been obtained without using many more parameters than are actually 
needed to obtain acceptably accurate results. In the problem under study, it was not pos- 
sible because of computer time limitations to increase the number of adjustable param- 
eters much beyond the point of apparent convergence. Thus, auxiliary tests must be used 
to check the results. 

Numerical results from the exact series solution given in table I can be compared 
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with the variational solution resvilts in table in for M values of 1, 2, 5, 10, and 20. The 
differences between these results are at most 1 part in the third place. To obtain a bet- 
ter feel for the magnitude of the error, the percentage difference between the exact and 
variational solutions is shown in figure 9 as a function of (R - l)y at M = 10. The differ- 
ences are determined for the average fluid velocity, induced potential difference, and 
center velocity. Figure 9 is based on results computed to 5 places which were later 
rounded to 3 places for presentation in tables I and HI. 

The figure shows that the center velocity has the largest error and the average veloc- 
ity has the least error. In fact, the error in the average velocity is lower by more than 
an order of magnitude from the error in the center velocity. This phenomenon was ex- 
pected since the stationary value of a functional, which in this case is proportional to the 
average velocity, is typically an order of magnitude more accurate than the error in the 
trial functions. The error in the induced potential difference also turns out to be 

about an order of magnitude less than the error in center velocity. Thus, the error in the 



Figure 9. - Error in solution from variational principle and Ritz technique 
at Hartmann number M = 10. 
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flowmeter sensitivity and distortion factor is about an order of magnitude less than the 
error in the center velocity. The same relations between the various errors were also 
fovmd for results at M values of 1, 2, and 5 although the magnitude of the errors de- 
creased as M decreased. 

Figure 9 also shows that the errors decrease as (R - l)y increases . Thus , the non- 
conductive pipe wall gives the largest error , and the thickest most highly conductive pipe 
wall gives the least error. This phenomenon was also expected since it took the largest 
number of parameters to obtain convergence for the case of a nonconductive pipe wall. 
The same number of parameters required for convergence at a given M value for the 
nonconductive wall case were used for all conductive wall cases at the same value for M. 
Consequently, the accuracy of the results should improve as the quantity (R - l)y in- 
creases. 

The error in the results for higher M values can be estimated by assuming that the 
effects observed for low M are also true for high M (i.e. , the errors in sensitivity and 
distortion factor are about an order of magnitude less than the error in center velocity 
and that a nonconductive pipe wall gives the largest error) . 

The center velocity for the nonconductive pipe wall case is the only quantity that can 
be easily computed exactly at any value of M. It is given by (eq. (21a) with p = 0) 



Figure 10 shows the percentage error in the variational solution for the dimension- 



Figure 10. - Error in solution for center velocity for a nonconductive pipe wali (solution obtained using vari- 
ational principle and Ritz technique). 
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less center velocity with a nonconductive pipe wall as a fxmction of M. As shown, the 
error increases with increasing M and is slightly less than 4 percent at M = 1000. The 
flattening effect in the error curve at M values of 20 to 100 is due to the large increase 
that was made in the number of adjustable parameters at these M values (see table II) 

If the error relations observed at low M are valid for hi^ M, the error in flowmeter 
sensitivity and distortion factor should be about 0. 4 percent at M = 1000 for a noncon- 
ductive pipe wall and decrease from this value as M decreases and (R - l)y increases. 

Other checks on the results can be made at high M. It is known from the general 
theory of magnetohydrodynamic channel flow that if the channel wall is conductive, the 
velocity profile is constant over the channel cross section almost up to the wall where it 
drops rapidly to zero (refs. 13 and 14). Thus, the distortion factor for the circular 
channel under consideration should approach 1 as M becomes large which is in agree- 
ment with the computed results. 

If the channel wall is nonconductive, the velocity profile remains distorted as M be- 
comes large. Although all the mathematical steps can not be formally justified, it has 
been shown by Shercliff (ref. 13) and Gold (ref. 9) that the flowmeter sensitivity and dis- 

o 

tortion factor are asymptotic to 37T /32 or 0.925 at high M. This value differs from the 
computed value of 0.924 at M = 1000 by only 0.16 percent. 

Another possible way of checking the results is to examine the energy balance. The 
rate at which energy is being supplied to the moving fluid by the pressure source should 
equal the rate at which energy is beii^ dissipated because of viscous losses in the fluid 
and ohmic losses in the fluid and wall. However, it can be shown that trial functions of 
the form given by equations (32) and (34) result in solutions which satisfy the energy bal- 
ance equation exactly so that this check is not an independent test. 

The approximate solutions can be tested, however, by determining how well they 
satisfy the basic equations (8a) to (8g). It is only necessary to examine equations (8a), 
(8b), and (8f) since the trial functions were selected to satisfy the remaining equations 
identically. The residuals, or differences between the left and right sides of these equa- 
tions were computed at selected points over the pipe cross section p ^ 1 for equations 
(8a) and (8b) and along the fluid - pipe wall interface p = 1 for equation (8f). 

For low values of M, it was found that if the residuals were small everywhere, then 
the approximate solution was almost exact as should be expected. Even if the residuals 
became large, however, which occasionally occurred near p = 1 for equation (8b), the 
approximate solution was still quite accurate. It was also observed that the residuals at 
all values of M became smaller as (R - l)y increased which indicates that the error at 
all values of M is reduced as (R - l)y is increased. 

Rather than attempt to determine a relation between the size of the residuals and the 
error in the solution, which is a very difficult task, it was decided to test the solutions 
as follows. Using the approximate solution for the velocity, the potential was computed 
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by solving equations (8a), (8c), (8e), (8f), and (8g) exactly. The potential difference 
for this solution was then compared with that obtained directly from the variational solu- 
tion. It was foimd that for M values of 20 or less the two solutions differed at most by 
0.01 percent. For M values of 50 to 1000 the two solutions differed typically by 0.05 
percent and at most by 0.15 percent. 

The values shown in table III for the potential flowmeter sensitivity, and 

distortion factor were actvially computed by solving for the potential exactly usir^ the ap- 
proximate variational solution for the velocity. The difference between the solutions ob- 
tained this way and those obtained directly from the variational solution, however, are at 
most 2 parts in the third place. 

The same procedure was also applied to equations (8b) and (8d) by solving them ex- 
actly for V usii^ the approximate solution for W. The difference between the average 
value of this solution for V and the average found from the stationary value of the func- 
tional was at most a few parts in the fifth place. 

As a final test, the method of solution using the Ritz technique was amended so that 
the trial functions for W exactly satisfied equations (8a), (8c), (8e), (8f), and (8g) for 
any given trial fvinction for V. The advantage of this method is that for any choice for the 
trial function for V, the trial function for W is constrained so as to satisfy all the basic 
equations identically except for equation (8b). Thus, all the adjustable parameters can be 
concentrated in the trial function for V. The sole purpose of the variational principle 
and Ritz technique then reduces to generating the best approximate solution for equa- 
tion (8b). This method of solution required a considerable amount of computer time so 
that it only was used for spot checking . The results obtained by this method differed very 
little from those given previously, however, indicating that the previous values are quite 
accurate. For example, at M = 200, R = 1.00, and y = 0 the previous solution for the 
center velocity which was in error by 2.2 percent was reduced to. 1.6 percent using this 
amended technique and 26 adjustable parameters compared with 40 parameters used pre- 
viously (see table II). The corresponding changes in average velocity, output potential, 
flowmeter sensitivity, and distortion factor were 0. 07, 0.16, 0.09, and 0.09 percent, 
respectively. 

On the basis of all these tests it is felt that the error in the main quantities of inter- 
est, namely the flowmeter sensitivity and distortion factor, are of the order of a few 
tenths of a percent at the higher Hartmann numbers which are commonly found in liquid 
metal applications. This error is comparable to the error in determining the pipe shimt- 
ing factor. Since the uncertainty in the electrical conductivity for most of the liquid 
metals is of the order of ±1 percent, the pipe shunting factor can only be determined 
typically to within ±0.25 percent. In addition, there are the errors associated with mea- 
suring the applied static magnetic field and output potential. However, these errors can 
often be held to less than ±0. 1 percent. 
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CONCLUDING REMARKS 


The theory of circular, transverse -field electromagnetic flowmeters has been ex- 
tended to include the effects of distortion from axial symmetry in the fluid velocity pro- 
file. A general formxilation was developed to relate the flowmeter sensitivity to any given 
but arbitrary velocity profile. It was found that flowmeters with a thick, highly conduc- 
tive wall are the least sensitive to uncertainties in the velocity profile. The analysis was 
then extended to determine the effect of velocity profile distortion by the flowmeter itself. 
The fully developed profile was computed usii^ a variational principle and the Ritz tech- 
nique. It was found that if the distortion from axial symmetry in the fully developed pro- 
file were not considered in computing the flowmeter sensitivity, then errors of up to 
8.2 percent would result in the flow measurement if the pipe wall were nonconductive . 

The error is reduced if the pipe wall is conductive, but it is still typically 2 to 5 percent 
for most liquid metal applications. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 28, 1971, 

120-27. 
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SYMBOLS 

expansion coefficient in series for f (eq. (20a)) 
expansion coefficient in trial function for f^ (eq. (34a)) 
pipe inner radius 
unit vectors 

dimensionless induced magnetic field (= B^/B^) 
magnetic field 

expansion coefficient in series for g (eq. (20b)) 
magnitude of B^ 

applied uniform static magnetic field 
induced axial magnetic field 
outer radius of pipe 

expansion coefficient for W in pipe wall (eq. (21c)) 
expansion coefficient in trial function for g^^ (eq. (34b)) 
electric field 

fvinctional for MHD flow (eq. (30)) 
flowmeter distortion factor (eq. (29)) 
stationary value of F 
auxiliary fvmction (eqs. (18a) and (18b)) 
radial fimction in expansion for V (eq. (32a)) 

Green's function 

radial fimction in expansion for G (eq. (12a)) 
auxiliary function (eqs. (18a) and (18b)) 
radial function in expansion for W (eq. (32b)) 

j-t- 

number of radial terms in n*"" angular term in expansion for V (eq. (34a)) 
electric current density 

number of radial t<"r ms in n “ angular term in expansion for W(eq. (34b)) 
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L 

M 
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P 

R 

m 

Re 

^®crit 

r,e,z 

S 

®nk 

Tn 

U 

^AD’^BC 

V 

V 
V„ 


m 


V, 


zm 

W 

Wad.Wbc 

x,y,z 

«ni 

^nj 


flowmeter weighting factor 
entry ler^h 
Hartmann number (= 
number of angular terms 


dimensionless axial pressure gradient 


- -a^ 3p\ 
V ^o 


pressure in fluid 

ratio of pipe outer -to-inner radii (= b/a) 
magnetic Reynolds number (= 

Reynolds number (= Piu^Vq/j/) 
critical Reynolds number 
cylindrical coordinates (fig. 2) 

flowmeter sensitivity (.U^/2aB^V^^=W^/2V„) 
infinite series (eq. (23f)) 
infinite series (eq. (23e)) 
induced potential 

induced potential difference between points A-D or 
dimensionless axial velocity (= V^/V^) 
fluid velocity 

dimensionless center velocity 
dimensionless average velocity 
characteristic fluid velocity 
axial fluid velocity 
average velocity 

dimensionless induced potential (= U/B^aV^) 

dimensionless potential difference between points A-] 

rectangular coordinates (fig. 2) 

exponent in expansion for f^ (eq. (34a)) 

exponent in expansion for g^ (eq. (34b)) 

ratio of wall to fluid conductivity (= 


-C (fig. 2) 


or B-C (fig. 2) 
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Kroniker delta function (=1 if i = j,=0 if i;^]) 
(= 2 if n = 0, = 1 if n 7 ^ 0) 
electric permittivity of free space 
17 fluid viscosity 

Pq magnetic permeability of free space 

p dimensionless radius (= r/a) 

Pg electric charge density 

mass density of fluid 

O’j electrical conductivity of fluid 

0'^ electrical conductivity of pipe wall 

^1 ’ ' ■ ■ ’ ^n adjustable prameters 
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APPENDIX B 


FLOWMETER PARAMETERS 

The purpose of this section is to acquaint the reader with some typical values for the 
flowmeter parameters that are encountered in liquid metal applications. First, consider 
the ratio of outer -to- inner pipe radii R. Figure 11 shows R as a fimction of the pipe 
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Figure 11. - Ratio of outer-to-inner radii for commercial steel pipes. 


inner radius a for schedule 20, 40, 60, and 80 commercial steel pipes. Despite the 
large variation in diameters and wall thicknesses for these pipes, the points on the 
curve R against a fall within a narrow band. Because of the large number of applications 
for pipes in the inner radius rai^e 0.01 to 0,1 meter, the R values used in the compu- 
tations were concentrated in the corresponding range of 1 . 00 to 1 . 35. In addition to the 
various schedule pipes, tubing is frequently used in liquid metal applications; particu- 
larly thin wall tubii^ in proposed space applications where weight is an important factor. 
The R values for this tubir^ are generally in the range 1.02 to 1.10. 

Next, consider the ratio of the pipe wall to fluid electrical conductivity y. Fig- 
iu*e 12 shows this ratio as a function of temperature for liquid mercury, lithium, sodium, 
potassium, and NaK for pipe walls made of 316 stainless steel and Till , a tantalum base 
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2 . 5 , 



Figure 12. - Ratio of wall to fluid electrical conductivity for common liquid metals and pipe materials. 


alloy. The conductivities of the liquid metals were obtained from references 23 and 24 
and the conductivities of the pipe wall materials were obtained from references 25 and 26. 
As shown in figure 12, the values of y encountered in practice are rather \mtformly dis- 
tributed over the range 0 to 2. 
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TABLE I. - RESULTS FROM EXACT SERIES SOLUTION 


(a) Hartmann number, M = 1 


Pipe outer -to- 
inner radii 
ratio, R 


1.00 

1.05 

1.15 

1.25 

1.35 


1.00 

1.05 

1.15 

1.25 

1.35 


Pipe wall to fluid conductivity ratio, y 


0.1 


0.5 


1.0 


1.5 


2.0 


Average velocity, V^^/P^ 


1.22x10 


-1 


1.22X10 


-1 


1.22x10 
1.22 
1.21 
1.21 


-1 


1.22x10 
1.21 
1.20 
1.19 


-1 


1.21X10 
1.20 
1.19 
1.18 


-1 


1.21x10 
1.19 
1.18 
1.17 


-1 


Induced potential, 


2.45X10 


-1 


2.43x10 
2.39 
2.33 
2.27 


-1 


2.38x10 
2.25 
2.13 
2.01 


-1 


2.32x10 
2.10 
1.92 
1.76 


-1 


2.26x10' 

1.96 

1.74 

1.57 


Sensitivity, S 


1.00 

0.998 




0.992 

0.973 

0.951 

1.05 



1.J5 

.975 

.925 

.868 


1.25 

.953 

.878 

.799 


1.35 

.928 

.834 

.740 





Distortion factor, 

1.00 

0.998 




0.998 

0.998 

0.998 

1.05 



1.15 

.998 

.999 

.999 


1.25 

.999. 

.999 

.999 


1.35 

.999 

.999 

.999 



0.929 

.819 

.733 

.665 


0.999 


2.20X10 
1.85 
1.60 
1.42 


-1 


0.909 

.774 

.677 

.604 


' 0.999 




TABLE I. - Continued. RESULTS FROM EXACT SERIES SOLUTION 


(c) Hartmann number, M = 5 


Pipe outer -to- 


Pipe wall to fluid conductivity ratio, y 


inner radii 
ratio, R 







0 

0.1 

0.5 

1.0 

1.5 

2.0 




Average velocity, V^/P 

o 


1.00 

8.96x10 ^ 






8.87x10"^ 

8.51x10"^ 

8.11x10"^ 

7.77x10"^ 

7.47X10”^ 

1.05 



1.15 

8.69 

7.82 

7.03 

6.46 

6.02 


1.25 

8.55 

7.33 

6.38 

5.76 

5.32 


1.35 

8.43 

6.97 

5.95 

5.32 

4.90 




Induced potential, 


1.00 

1. 75x10"^ 






1.72x10"^ 

1.62x10"^ 

1.51X10"^ 

1.42X10"^ 

1.33x10”^ 

1.05 



1.15 

1.66 

1.42 

1.20 

1.04 

9.19x10"^ 


1.25 

1.60 

1.27 

1.01 

8.33x10"^ 

7.12 


1.35 

1.54 

. 

1.15 

8.70x10"^ 

7.00 

5.86 





Sensitivity, S 



1.00 

0.975 






0.969 

0.952 

0.931 

0.911 

0.892 

1.05 



1.15 

.955 

.907 

.854 

.806 

.763 


l.,25 

.935 

.864 

.788 

.724 

.669 


1.35 

.913 

.823 

.731 

.658 

.598 





Distortion factor, 



1.00 

0.975 






0.975 

0.976 

0.978 

0.979 

0.980 

1.05 



1.15 

.977 

.980 

.982 

.983 

.985 


1.25 

.980 

.983 

.985 

.986 

.987 


1.35 

Q Q9 

.985 

.987 

.988 

.989 
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TABLE I. - Continued. RESULTS FROM EXACT SERIES SOLUTION 


(d) Hartmann number, M = 10 


Pipe outer -to- 

Pipe wall to fluid conductivity ratio, y 

inner radii 







ratio, R 

0 

0.1 

0.5 

1.0 

1.5 

2.0 



Average velocity, V^/P 

0 


1 00 

6.01x10"^ 






1.05 

5. 84x10'^ 

5.25x10"^ 

4.68x10"^ 

4. 24x10"^ 

3.90x10"^ 

1.15 


5.55 

4.30 

3.45 

2.93 

2.59 

1.25 


5.31 

3.75 

2.87 

2.40 

2.11 

1.35 


5.12 

3.39 

2.54 

2.11 

1.86 



Induced potential, W^j^/P^ 


1 00 

1.15x10'^ 






1.05 

1. 11X10"^ 

9.80x10‘2 

8. 57X10"^ 

7.61X10'^ 

6.85X10“^ 

1.15 


1.04 

7.69 

5.82 

4.68 

3.92 

1.25 


9.77X10'^ 

6.40 

4.48 

3.45 

2.80 

1.35 


9.23 

5.52 

3.68 

2.76 

2.21 




Sensitivity, S 



1 . 00 

0.954 






1.05 

0.949 

0.934 

0.915 

0.897 

0.879 

1.15 


.936 

.894 

.844 

.798 

.756 

1.25 


.920 

.853 

.781 

.718 

.665 

1.35 


.900 

.814 

.726 

.654 

.595 


Distortion factor, Fj^ 

1 , 00 

0.954 






1.05 

0.955 

0.958 

0.961 

0.964 

0.966 

1.15 


.959 

.965 

.970 

.974 

.976 

1.25 


.964 

.971 

.976 

.979 

.981 

1.35 


.968 

.975 

.980 

.982 

.984 


TABLE I. - Concluded. RESULTS FROM EXACT SERIES SOLUTION 


(e) Hartmann number, M = 20 


Pipe outer -to- 


Pipe wall to fluid conductivity ratio, y 


inner radii 






1 

ratio, R 

0 

0.1 

0.5 

1.0 

1.5 

2.0 




Average velocity, V^/P^ 


1 00 

3.53X10"^ 






1.05 

3.30X10"^ 

2.62x10"^ 

2.11X10"^ 

1.78X10'^ 

1.55x10"^ 

1.15 


2.95 

1.82 

1.28 

1.02 

8.59X10"^ 

1.25 


2.69 

1.46 

9.87X10”^ 

7.78x10"^ 

6.60 

1.35 


2.51 

1.25 

8.36 

6.61 

5.65 



Induced potential, 


1.00 

6.63X10"^ 











1.05 


6.17X10"^ 

4.85X10"^ 

3.84X10"^ 

3.18X10“^ 

2.71X10”^ 

1.15 


5.45 

3.24 

2.16 

1.62 

1.30 

1.25 


4.91 

2.48 

1.54 

1.12 

8.79x10"^ 

1.35 


4.48 

2.03 

1.21 

8.66X10"^ 

6.74 




Sensitivity, S 



1 . 00 

0.938 












1.05 


0.935 

0.924 

0.909 

0.892 

0.876 

1.15 


.925 

.889 

.842 

.798 

.757 

1.25 


.911 

.851 

.781 

.719 

.666 

1.35 


.893 

.813 

.726 

.655 

.596 




Distortion factor, 



1.00 

0.938 











1.05 

— 

0.941 

0.948 

0.954 

0.959 

0.963 

1.15 


.947 

.960 

.969 

.974 

.977 

1.25 


.955 

.968 

.976 

.980 

.982 

1.35 


.961 

.974 

.981 

.984 

.986 
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TABLE n. - NUMBER OF ADJUSTABLE PARAMETERS IN TRIAL 
FUNCTIONS FOR VELOCITY AND POTENTIAL 


Hartmann 

Number of 

Number of radial terms in each angular term 

number , 
M 

angular 
terms , 


Velocity 



Potential 



N 


h 

^3 

I4 


^2 

J3 

J4 

1 

2 

3 

2 

- 

- 

3 

2 

~ 

- 

2 

2 

3 

2 

- 

- 

3 

2 

- 

- 

5 

2 

3 

2 

- 

-- 

3 

2 

“ 

- 

10 

2 

3 

2 

- 

- 

3 

2 

- 

- 

20 

2 

3 

2 

- 

- 

3 

2 

- 

- 

50 

4 

4 

2 

3 

3 

3 

2 

2 

- 

100 

4 

4 

3 

3 

3 

3 

2 

2 

- 

200 

4 

4 

3 

3 

3 

3 

2 

2 

- 

500 

4 

4 

3 

3 

3 

3 

2 

2 

- 

1000 

4 

4 

3 

3 

3 

3 

2 

2 

- 
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TABLE ni. - RESULTS FROM VARIATIONAL PRINCIPLE AND RITZ TECHNIQUE 


(a) Hartmann number, M = 1 


Pipe outer -to- 

Pipe wall to fluid conductivity ratio, y 

inner radii 









ratio, R 

0 

0.1 

0.5 

1.0 


1.5 


2.0 


Average velocity, V^/P^ 

1 00 

1 . 22 X 10 "^ 








1.05 

1 . 22 X 10 "^ 

1 . 22 x 10 "^ 

1 . 22 X 10 "^ 

1 . 

21 X 10 ’^ 

1 . 

21 x 10 "^ 

1.15 


1.22 

1.22 

1.21 

1 . 

20 

1 . 

19 

1.25 


1.22 

1.21 

1.20 

1 . 

19 

1 . 

18 

1.35 


1.22 

1.21 

1.19 

1 . 

18 

1 . 

17 


Induced potential, 

1 00 

2 . 45 X 10 "^ 




1 


1.05 

2 . 43 X 10 "^ 

2 . 38 X 10 "^ 

2 . 32 X 10 "^ 

2 . 

26 X 10 "^ 

2 . 20 X 10 '^ 

1.15 


2.39 

2.25 

2.10 

1 . 

96 

1 . 

85 

1.25 


2.33 

2.13 

1.92 

1 . 

74 

1.60 

1.35 


2.27 

2.01 

1.76 

1 . 

57 

1.42 


Sensitivity, S 

1 . 00 

0.998 




■ 


1.05 

0.992 

0.973 

0.951 

0.929 

0.909 

1.15 


.975 

.925 

.868 

.819 

.774 

1.25 


.953 

.878 

.799 

.733 

.677 

1.35 


.928 

.834 

.740 

.665 

.604 


Distortion factor, 

1 nn 

0.998 




1 

1 









1.05 


0.998 

0.998 

0.998 

0.999 

0.999 

1.15 


.998 

.999 

.999 





1.25 


.999 

.999 

.999 





1.35 


.999 

.999 

.999 






62 



TABLE m. - Continued. RESULTS FROM VARIATIONAL PRINCIPLE AND 


RITZ TECHNIQUE 


(b) Hartmann number, M = 2 


Pipe outer -to- 


Pipe wall to fluid conductivity ratio, y 


inner radii 

1 






ratio, R 

0 

0.1 

0. 5 

1.0 

1.5 

2.0 



Average velocity, V^/P 

0 


1.00 













1.05 


1. 16X10"^ 

1.15x10“^ 

1.13X10"^ 

1.12x10”^ 

1. 11X10"^ 

1.15 


1.15 

1.13 

1.10 

1.07 

1.05 

1.25 


1.15 

1.11 

1.07 

1.04 

1.02 

1.35 


1.14 

1.09 

1.05 

1.02 

9.90X10”^ 



Induced potential, 


1.00 

2.30x10"^ 











1.05 


2.28X10"^ 

2.22X10"^ 

2.15X10"^ 

2.08x10"^ 

2.02x10"^ 

1.15 


2.24 

2.07 

1.90 

1.75 

1.63 

1.25 


2.18 

1.94 

1.70 

1.52 

1.37 

1.35 


2.12 

1.82 

1.55 

1.35 

1.19 




Sensitivity, S 



1.00 

0.994 











l.,05 


0.988 

0.969 

0.947 

0.926 

0.905 

1.15 


.971 

.921 

.866 

.816 

.772 

1.25 


.950 

.875 

.797 

.731 

.676 

1.35 


.925 

.832 

.738 

.664 

.603 




Distortion factor, Fj^ 



1 . 00 

0.994 






1,05 

0.994 

0.994 

0.994 

0.995 

0.995 

1.15 


.994 

.995 

.995 

.996 

.996 

1.25 


.995 

.996 

.996 

.996 

.997 

1.35 


.996 

.996 

.997 

.997 

.997 



II 


I iiiiiii I iiiiii II I nil I 


■■■■I I II 
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TABLE m. - Continued. RESULTS FROM VARIATIONAL PRINCIPLE AND 

RITZ TECHNIQUE 


(c) Hartmann number, M = 5 


Pipe outer -to- 
inner radii 
ratio, R 

Pipe wall to fluid conductivity ratio, y 

0 

0.1 

0.5 

1.0 

1.5 

2.0 




Average velocity, 


1.00 

8.96X10"^ 






8.86X10”^ 

8.51X10"^ 

8.11X10"^ 

7.77X10"^ 

7.47X10'^ 

1.05 



1.15 

8.69 

7.82 

7.03 

6.46 

6.02 


1.25 

8.56 

7.33 

6.38 

5.76 

5.32 


1.35 

8.42 

6.97 

5.95 

5.32 

4.90 




Induced potential, 


1.00 

1.75x10"^ 






1.72x10"^ 

1.62X10"^ 

1. 51X10’^ 

1.42X10"^ 

1.33X10"^ 

1.05 



1.15 

1.66 

1.42 

1.20 

1.04 

9.19X10"^ 


1.25 

1.60 

1.27 

1.01 

8.33x10‘2 

7.12 


1.35 

1.54 

1.15 

8.69X10'^ 

7.00 

5. 86 





Sensitivity, S 



1.00 

0.974 





1 

0.969 

0.952 

0.931 

0.911 

0. 891 

1.05 



1.15 

.954 

.907 

.854 

.806 

.763 


1.25 

.935 

.864 

.788 

.724 

.669 


1.35 

.913 

.822 

.731 

.658 

.598 





Distortion factor , 



1.00 

0.974 






0.975 

0.976 

0.977 

0.978 

0.979 

1.05 



1.15 

.977 

.979 

.982 

.983 

.984 


1.25 

.980 

.982 

.985 

.986 

.987 


1.35 

.982 

.985 

.987 

.988 

.989 
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TABLE m. - Continued. RESULTS FROM VARIATIONAL PRINCIPLE AND 


RITZ TECHNIQUE 


(d) Hartmann number, M = 10 


Pipe outer -to- 


Pipe wall to fluid conductivity ratio, y 


inner radii 







ratio, R 

0 

0.1 

0.5 

1.0 

1.5 

2.0 




Average velocity, V^^/P^ 


1 00 

6.01X10“^ 






1.05 

5. 83X10"^ 

5.24x10'^ 

4.68X10“^ 

4.24X10'^ 

3.90X10"^ 

1.15 


5.54 

4.30 

3.45 

2.93 

2.59 

1.25 


5.31 

3.75 

2.87 

2.40 

2.11 

1.35 


5.12 

3.39 

2.54 

2.11 

1.86 



Induced potential, 


1 00 

1.14x10'^ 






1.05 

1. 11X10"^ 

9.79X10"^ 

8. 56X10"^ 

7.61X10"^ 

6.85X10"^ 

1.15 


1.04 

7.69 

5.82 

4.68 

3.92 

1.25 


9.77X10'^ 

6.40 

4.48 

3.45 

2.80 

1.35 


9.22 

5.52 

3.68 

2.76 

2.21 




Sensitivity, S 



1 . 00 

0.952 






1.05 

0.948 

0.933 

0.914 

0.896 

0.878 

1.15 


.936 

.893 

.843 

.798 

.756 

1.25 


.919 

.853 

.780 

.718 

.665 

1.35 


.900 

.814 

.726 

.654 

.595 




Distortion factor, Fj^ 



1.00 

0.952 











1.05 


0.954 

0.957 

0.960 

0.963 

0.965 

1.15 


.958 

.964 

.970 

.973 

.976 

1.25 


.963 

.970 

.976 

.979 

.981 

1.35 


.968 

.975 

.980 

.982 

.984 


TABLE m. - Continued. RESULTS FROM VARIATIONAL PRINCIPLE AND 


RITZ TECHNIQUE 


(e) Hartmann number, M = 20 


Pipe outer -to- 


Pipe wall to fluid conductivity ratio, y 


inner radii 
ratio, R 

0 

0.1 

P 0.5 

1.0 

1.5 

2.0 




Average velocity , V^/P^ 


1.00 

3.53X10"^ 






3.30X10"^ 

2.62X10"^ 

2.11X10"^ 

1.78X10"^ 

1.55X10"^ 

1.05 



1.15 

2.94 

1.82 

1.28 

1.02 

8.59X10"^ 


1.25 

2.69 

1.46 

9.87X10"^ 

7.78X10"^ 

6.60 


1.35 

2.51 

1.25 

8.36 

6.61 

5.65 




Induced potential, 


1.00 

6.62x10"^ 






6.17X10"^ 

4.85x10"^ 

3.84X10"^ 

3.18x10"^ 

2.71X10'^ 

1.05 



1.15 

5.45 

3.24 

2.16 

1.62 

1.30 


1.25 


4.91 

2.48 

1.54 

1.12 

Q r70s/1 


o. lyAiu 
6.74 

1.35 

4.48 

2.03 

1.21 

8.66X10"^ 





Sensitivity, S 



1.00 

0.937 






0.935 

0.924 

0.909 

0.893 

0. 876 

1.05 



1.15 

.925 

.890 

.843 

.799 

.758 


1.25 

.911 

.852 

.781 

.720 

.666 


1.35 

.893 

.814 

.727 

.655 

.596 





Distortion factor, 



1.00 

0.937 






0.940 

0.948 

0.954 

0.959 

0.963 

1.05 



1.15 

.947 

.961 

.969 

.974 

.978 


1.25 

.955 

.969 

.976 

.980 

.983 


1.35 

.961 

.974 

.981 

.984 

.986 
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TABLE ni. - Continued. RESULTS FROM VARIATIONAL PRINCIPLE AND 


RITZ TECHNIQUE 


(f) Hartmann number, M = 50 


Pipe outer -to- 


Pipe wall to fluid conductivity ratio, y 


inner radii 

— 



— 



ratio, R 

0 

1 0.1 

0.5 

1.0 

1.5 

2.0 




Average velocity, V^/P^ 


1.00 

I . 57 X 10 “^ 






1.05 

1 . 31 X 10 '^ 

8 . 01 X 10 "^ 

5 . 48 X 10 ’^ 

4 . 22 x 10 ’^ 

3 . 46 X 10 ’^ 

1.15 


1.01 

4.37 

2.69 

2.01 

1.64 

1.25 


8 . 41 X 10 "^ 

3.18 

1.93 

1.45 

1.20 

1.35 


7.34 

2.60 

1.58 

1.21 

1.01 



Induced potential, 


1.00 

2 . 91 X 10 "^ 

' ■ ‘ 





1.05 

2 . 43 X 10 “^ 

1 . 48 x 10 

I.OOXIO‘2 

7 . 60 X 10 '^ 

6 . 12 X 10 "^ 

1.15 


1.87 

7 . 84 X 10 "^ 

4 . 57 X 10 "^ 

3.23 

2.50 

1.25 


1.54 

5.47 

3.04 

2.11 

1.61 

1.35 


1.32 

4.26 

2.31 

1.59 

1.21 




Sensitivity, S 



1.00 

0.927 






1.05 

0.928 

0.926 

0.915 

0.901 

0.885 

1.15 


.924 

.897 

.851 

.805 

.764 

1.25 


.913 

.859 

.787 

.725 

.671 

1.35 


.896 

.820 

.732 

.659 

.599 




Distortion factor , Fj^ 



1.00 

0.927 






1.05 

0.934 

0.950 

0.961 

0.968 

0.972 

1.15 


.946 

.969 

.978 

.983 

.985 

1.25 


.956 

.977 

.984 

.987 

.989 

1.35 


.964 

.982 

.988 

.990 

.991 



TABLE m. - Continued. RESULTS FROM VARIATIONAL PRINCIPLE AND 


RITZ TECHNIQUE 


(g) Hartmann number, M = 100 


Pipe outer -to- 


Pipe wall to fluid conductivity ratio, y 


inner radii 







ratio, R 

0 

0.1 

0.5 

1.0 

1.5 

2.0 




Average velocity, V^/P^ 


1 00 

8.15X10"^ 






1.05 

5.78x10“^ 

2.73X10'^ 

1.68X10"^ 

1.23X10'^ 

9.81X10"^ 

1.15 


3.80 

1.28 

7. 40X10"^ 

5.39X10'^ 

4.34 

1.25 


2.92 

8.93X10’'^ 

5.18 

3.83 

3.14 

1.35 


2.43 

7.13 

4.19 

3.15 

2.62 



Induced potential, 


1.00 

l.SlxlO'^ 











1.05 


1.08X10'^ 

5.12X10"^ 

3.11X10"^ 

2.24X10'^ 

1.75X10'^ 

1.15 


7.08X10'^ 

2.33 

1.27 

8.74X10’'^ 

6.66X10"^ 

1.25 


5.38 

1.55 

8.20x10"'^ 

5.58 

4.23 

1.35 


4.40 

1.18 

6.15 

4.17 

3.15 




Sensitivity, S 



1 . 00 

0.925 






1.05 

0.932 

0.936 

0.926 

0.910 

0.893 

1.15 


.932 

.906 

. 857 

.811 

.768 

1.25 


.921 

.866 

.792 

.728 

.673 

1.35 


.904 

.825 

.735 

.661 

.601 




Distortion factor, 



1 . 00 

0.925 






1.05 

0.937 

0.960 

0.972 

0.978 

0.982 

1.15 


.954 

.979 

.986 

.989 

.991 

1.25 


.965 

.985 

.990 

.992 

.993 

1.35 


.972 

.989 

.992 

.994 

.995 
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TABLE m. - Continued. RESULTS FROM VARIATIONAL PRINCIPLE AND 


RITZ TECHNIQUE 


(h) Hartmann number, M = 200 


Pipe outer -to- 


Pipe wall to fluid conductivity ratio, y 


inner radii 







ratio, R 

0 

0.1 

0.5 

1.0 

1.5 

2.0 




Average velocity, V^/P 

o 


1.00 

4.15x10"^ 






1.05 

2.27X10"^ 

8.29X10''^ 

4.72X10’^ 

3.35X10"^ 

2.62x10"'^ 

1.15 


1.25 

3.50 

1.95 

1.40 

1.12 

1.25 


8.99X10“^ 

2.37 

1.34 

9.83X10"® 

8.02x10"® 

1.35 


7.21 

1.87 

1.08 

8.05 

6.67 



Induced potential, W^^/P^ 


1.00 

7.67X10’^ 






1.05 

4.26X10"^ 

1.57X10"^ 

8. 84X10"^ 

6.16X10"^ 

4.72X10"“* 

1.15 


2.36 

6.41X10"^ 

3.36 

2.28 

1.72 

1.25 


1.67 

4.14 

2.13 

1.44 

1.08 

1.35 


1.31 

3.10 

1.59 

1.07 

8.04X10"® 




Sensitivity, S 



1 . 00 

0.924 






1.05 

0.939 

0.948 

0.936 

0.919 

0.900 

1.15 


.942 

.914 

.862 

.814 

.771 

1.25 


.930 

.872 

.795 

.730 

.675 

1.35 


.912 

.829 

.737 

.663 

.602 




Distortion factor, Fj^ 



1.00 

0.924 






1.05 

0.945 

0.972 

0.983 

0.987 

0.989 

1.15 


.964 

.987 

.992 

.993 

.994 

1.25 


.975 

.992 

.994 

.995 

.996 

1.35 


.981 

.993 

.995 

.996 

.997 


TABLE III. - Continued. RESULTS FROM VARIATIONAL PRINCIPLE AND 


RITZ TECHNIQUE 
(i) Hartmann number, M = 500 


Pipe outer-to- 
inner radii 
ratio, R 


Pipe wall to fluid conductivity ratio, y 


0 

0.1 

0. 5 

1.0 

1.5 

2.0 




Average velocity, V^/P^ 


1.00 

1.68X10'^ 






5.45x10'^ 

1.52X10'^ 

8.15X10*® 

5.66X10*® 

4. 38X10"® 

1.05 



1.15 

2.46 

5.93X10"^ 

3.21 

2.29 

1.82 


1.25 

1.66 

3.94 

2.19 

1.60 

1.30 


1.35 

1.30 

3.08 

1.75 

1.30 

1.08 




Induced potential, 

P 

0 


1.00 

3.10x10'^ 






1.04X10'^ 

2.93X10'^ 

1.54X10*^ 

1.05X10*^ 

7.93X10*® 

1.05 



1.15 

4.71x10’'^ 

1.09 

5.57X10*® 

3.74X10*® 

2.81 


1.25 

3.13 

6.91x10*® 

3.50 

2.34 

1.76 


1.35 

2.39 

5.13 

2.59 

1.73 

1.30 





Sensitivity, S 



1.00 

0.924 






0.952 

0.965 

0.946 

0.926 

0.906 

1.05 



1.15 

.955 

.921 

.866 

.817 

.773 


1.25 

.941 

.876 

.798 

.732 

.677 


1,35 

Q 0/1 

.833 

.739 

.664 

.603 


. 




Distortion factor , 



1.00 

0.924 






0.958 

0.989 

0.994 

0.995 

0.995 

1.05 



1.15 

.978 

.995 

.996 

.997 

.997 


1.25 

.986 

.996 

.997 

.998 

.998 


1.35 

.994 

.997 

.998 

.998 

.998 
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TABLE m. - Concluded. RESULTS FROM VARIATIONAL PRINCIPT.E AND 


UlTZ TECHNIQUE 


(j) Hartmann number, M - 1000 


Pipe outer-to- 


Pipe wall to fluid conductivity ratio, y 


inner radii 
ratio, R 

0 

0.1 

0.5 

1.0 

1.5 

2.0 




Average velocity, P 

o 


1.00 

8.43X10“^ 






1. 64X10“'^ 

3.99x10“^ 

2.09x10'® 

1.44X10"® 

1.11X10"® 

1.05 



1.15 

6.67X10“^ 

1.51 

8.12x10"® 

5.76X10"® 

4. 57X10"® 


1.25 

4.39 

9.98X10"® 

5.52 

4.02 

3.26 


1.35 

3.38 

7.78 

4.41 

3.27 

2.71 




Induced potential, 1 



1.00 

1.05 

1. 56X10“^ 






3.15X10"'* 

7.75X10"® 

3.98X10"® 

2.67X10"® 

2.01X10"® 


1.15 

1,29 

2.79 

1.41 

9.42x10"® 

7.08x10"® 


1.25 

8.33X10"^ 

1.75 

8.81X10"® 

5.89 

4.42 


1.35 

6.27 

1.30 

6.52 

4.35 

3.27 



Sensitivity, S 

1.00 

0.924 






0.963 

0.972 

0.949 

0.928 

0.908 

1.05 



1.15 

.963 

.924 

.868 

.818 

.774 


1.25 

.948 

,877 

.799 

.733 

.677 


1.35 

.928 

.834 

.740 

.665 

.604 





Distortion factor , 



1.00 

0.924 






0.969 

0.997 

0.997 

0.997 

0.997 

1.05 



1.15 

.986 

.997 

.998 

.998 

,998 


1.25 

.993 

.998 

.998 

.999 

.999 


1,35 

.998 

.998 

.999 

.999 

.999 
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